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ABSTRACT 


The  electrical  resistivities  of  single  crystal  NaCl,  K1 ,  and  Csl  have 
been  measured  at  several  shock  pressures  in  the  rang"  from  120  to  300  kbai  . 
The  resistivity  of  NaCl  remained  >  9  x  103  fl-cm  at  225  kbar.  One  shot  in 
which  the  specimen  was  preheated  to  480°C  and  shocked  to  244  kbar  gave  an 
averaged  resistivity  of  470  f2-cm.  The  resistivity  of  KI  decreased  to  a 
value  between  2  and  9  O-cm  at  241  kbar;  there  is  considerable  scatter  in 
the  data  below  190  kbar.  Resistivity  measurements  on  Csl  were  performed 
parallel  and  perpendicular  lo  the  direction  of  shock  propagation;  the 
resistivity  decreased  to  2  0-  cm  at  277  kbar  in  the  first  case  and  to 
~0.2H-cm  in  the  latter.  The  data  suggest  a  dependence  of  the  conductivity 
on  crystallographic  direction.  The  data  are  consistent  neither  with 
electronic  conduction  caused  by  a  decrease  in  the  band  gap  dje  to  shock 
compression  of  the  lattice  nor  with  ionic  conduction  caused  by  melting  which 
may  occur  behind  the  shock  front. 
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I  INTRODUCTION 


In  the  past  few  years,  many  1 aboratori es  have  studied  the  electrical 
effects  associated  with  the  propagation  of  shock  waves  in  solids.  Included 
in  these  effects  are  the  depolarization  of  ferroelectric  ceramics,  the  con¬ 
ductivity  of  sulfur  and  other  materials,  and  the  polarization  of  certain 
materials.  The  motivation  for  some  of  this  work  has  been  the  search  for 
explosive- sensi tive  and  pressure-electric  transducers. 

Basic  study  of  some  of  these  phenomena  should  produce  data  and  insight 
concerning  the  fundamentals  of  the  shock  process  itself  and  should  also  pro¬ 
vide  tools  for  study  of  the  changes  induced  in  solids  by  shock  waves.  It 
is  obvious,  however,  that  a  project  of  modest  size  must  concentrate  its 
efforts  on  a  few  specific  areas  of  interest  if  it  is  not  to  be  completely 
exploratory  in  nature.  There  is  indeed  considerable  incentive  for  explor¬ 
atory  research — the  need  for  diagnostic  tools  for  pressi re  and  temperature 
measurement,  the  need  to  avoid  spurious  electrical  effects  in  present  elec¬ 
trical  measurements  utilizing  circuit  elements  which  are  shocked  in  the 
course  of  the  measurement,  and,  of  course,  scientific  curiosity. 

In  1956,  Alder  and  Christian  reported  some  semiquantitative  data  on 
several  ionic  solids  which  indicated  that  their  electrical  conductivity 
dramatically  increased  when  the  solids  were  compressed  by  strong  shock 
waves.1  A  study  of  this  effect  was  begun  two  years  ago  (under  BRL 
sponsorship)*  at  Stanford  Research  Institute2  for  the  following  reasons: 

1.  The  effect  was  not  predictable  based  on  present  knowledge. 

2.  Efforts  by  Griggs,  et  al .  3  and  Drickamer4  to  detect  the  effect 

in  static  high-pressure  apparatus  were  unsuccessful.  The  need  for  corrob¬ 
oration  of  the  effect  was  heightened  Hy  the  then- primi ti ve  state  of 
electrical  conductivity  measurements  in  shock  experiments. 

3.  If  the  effect  were  corroborated,  it  would  be  important  to  make 
the  measurements  more  quantitative  and  to  measure  the  pressure  and  tem¬ 
perature  dependence  of  the  conductivity  so  that  the  mechanism  could  be 
identi f ied. 

*  A  report  covering  the  first  year’s  work  will  be  referred  to  herein  as  Report  I.  The  contents  of  Report  I 
will  not  be  repeated  here  unless  necessary  for  completeness. 
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4.  The  interpretation  of  the  effect  by  Alder  and  Christian1  as  a 
metallic  transition  was  not  supported  by  the  theoretical  work  of  Flower 
and  March. 5  Alternative  interpretations  needed  to  be  examined. 

During  the  period  covered  by  Report  I,  the  effect  was  corroborated, 
the  measurements  were  made  more  quantitative,  and  a  greater  range  of  pres¬ 
sure  and  temperature  was  investigated.  An  alternative  interpretation  in 
terms  of  an  ionic  mechanism  was  suggested,  but  the  mechanism  was  not 
positively  identified. 
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II  MATERIALS 


The  materials  studied  were  NaCl,  KI ,  and  Csl ;  the  major  emphasis 
during  1964  has  been  on  Csl. 

The  specimens  of  NaCl  and  KI  were  prepared  from  optical  grade  crystals 
obtained  from  Harshaw  Chemical  Co.,  Cleveland,  Ohio.  The  NaCi  specimens 
were  single  crystals.  The  KI  specimens  were  poor  single  crystals  in  that 
they  contained  domains  of  slight  misorientation  (±1°  according  to  Harshaw). 
They  were,  however,  easily  cleavable.  The  KI  and  NaCl  were  cleaved  to  the 
desired  thicknesses,  which  for  the  present  experiments  varied  from  1  to 
2-1/2  mm.  Surface  irregularities  produced  by  cleaving  were  removed  by 
lapping,  so  that  final  parallelism  was  within  0.001  inch.  The  KI  crystals 
were  kept  desiccated  until  fired. 

Csl  crystals  were  obtained  from  Harshaw  and  from  Isotopes  Inc.  The 
specimens  from  Harshaw,  as  noted  in  Report  I,  were  polycrystalline  with 
apparent  domain  size  of  the  order  of  1  mm;  however,  the  crystals  used  in 
Shots  10,820  and  10,943  were  single  crystals  obtained  from  Harshaw.  The 
Isotope  crystals  were  single  crystal. 

Csl  crystals,  not  being  cleavable,  originally  were  cut  either  with  a 
diamond  saw  or  with  a  thin  blade  held  stationary  as  the  crystal  rotated  in 
a  lathe.  During  the  latter  part  of  the  work  the  Csl  crystals  were  cut  with 
a  wire  saw,  and  the  surfaces  lapped  flat  with  600  grit  paper  and  alcohol. 

All  specimens  were  circular  disks  25  mm  in  diameter  except  those  used 
in  the  transverse  geometry  shots. 
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Ill  EXPERIMENTAL  TECHNIQUE 


A.  SHOT  CONFIGURATION 

The  “longitudinal  geometry”  in  which  the  direction  of  the  current  is 
parallel  to  the  propagation  direction  of  the  shock  is  shown  in  Figure  1. 

The  driver  plate,  which  also  served  as  one  electrode,  was  2024T  aluminum 
for  all  shots.  The  backing  electrode  was  aluminum  for  Csl  specimens  and 
magnesium  for  NaCl  and  KI  specimens,  in  order  to  match  approximately  the 
shock  impedances  of  specimen  and  electrode.  This  electrode  was  usually 
22  mm  in  diameter;  thus  when  no  insulation  was  used  the  outer  1-1/2  mm  of 
the  upper  surface  of  the  crystal  and  the  curved  edge  were  free.  For  some 
shots,  the  sample  was  insulated  laterally  with  Dow  Corning  silicone  fluid 
No.  200  (electrical  grade — 200  cs). 

The  electrical  leads  were  connected  to  the  backing  electrode  at  least 
12  mm  from  the  specimen,  so  that  the  shock  would  not  reach  the  junction 
during  the  time  of  measurement  1  /isec  after  shock  transit).  For  the 
same  reason,  the  ground  leads  were  connected  at  the  top  of  a  30-mm-high 
aluminum  cylinder  which  rested  on  the  driver  plate.  This  also  avoided  a 
possible  thermal  emf  which  could  be  generated  at  the  junction  of  the  steel 
screw  with  the  aluminum  driver  (see  Report  I,  Part  I)  if  the  screw  were 
used  as  a  portion  of  the  ground  connection. 

On  each  shot  we  included  a  pressure  transducer  of  the  resi  st ance  wi  re 
type  (Manganin  wire  in  C-7  epoxy).6  This  served  not  only  to  monitor  the 
shock  strength  but  also  to  provide  an  excellent  time  reference  signal.  The 
transducer  signal  and  the  specimen  signals  were  displayed  on  the  two  beams 
of  a  Tektronix  551  oscilloscope  for  direct  comparison.  The  gauge  wire  was 
1  mm  from  the  epoxy-driver  plate  interface,  but  shock  arrival  at  the  inter¬ 
face  was  indicated  by  a  small  polarization  signal7  and  hence  no  correction 
for  shock  transit  in  the  C-7  was  necessary.* 


*  Some  of  the  shots  included  s  gauge  of  slightly  different  design  from  that  developed  snd  tested  exten¬ 
sively  ir  these  lsborstories.  The  sltered  design  wss  intended  to  give  grester  time  coverage— which  it 
did.  However,  the  indicated  peak  pressures  were  erratic  and  have  been  discounted.  The  records  from 
these  gauges  were  used  only  as  s  time  reference. 
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\  lew  slid  t  s  were  I  i  r ed  in  vvliiali  I  wn  s  |  >  •  *  <  I  tsia  ■  li  s  »i-i  c  >  I  ad  I  <<l  s  i  mil  I  - 
tuncouslx.  Iliese  experiments  i  it  wlii:h  one  spei  i  <n  <  - 1  ■  a  semL  I  y  ,s<  it<d 

.1  s  a  rout  rot  lor  .1  parameter  change  in  tin-  other. 

Pop  tin*  "  t  ransuT.sc  geometry”  (Figure  2),  t  li<*  <iirr«nt  ill  rci  linn  is 
perpendicular  to  the  direction  of  slioek  propagation.  The  sample  was  a 
nearly  rectangular  block,  6.35  x  4  x  25.4  mm,  which  was  supported  on  all 
sides  (but  not  the  ends)  by  the  same  material  as  that  of  the  sample. 

Teflon  sheet,  0.003  inch  thick,  was  used  as  electrical  insulation.  The 
electrodes  were  strips  of  0.001-inch  aluminum  foil,  6.35  mm  wide,  placed 
equidistant  from  the  ends  of  the  sample  (see  Figure  2).  A  brass  bar  was 
fastened  with  epoxy  to  each  of  the  s i de  -  suppo r t i ng  blocks  and  all  electrical 
connections  were  matte  to  these  Lars.  The  constant- current  measuring  cir¬ 
cuit  was  used  (see  Section  B  below).  The  sample,  supporting  blocks, 
electrodes,  and  Teflon  insulation  were  ail  assembled  and  held  tightly  with 
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a  clamp  while  the  edges  were  cemented  with  epoxy.  this  sample  assembly  was 
insulated  electrically  from  the  aluminum  driver  by  a  1  / 1- i  rich- th  i  ok  plate 
of  Hi-D  glass  (density  6.2  g/cm3).*  A  Siangan  i  n  -  wi  re  pressure  transducer  of 
the  kind  described  earlier  was  mounted  on  a  similar  glass  plate. 

The  transverse  geometry  technique  allows  a  measurement  of  the  sample 
resistance  as  a  function  of  time  during  shock  passage  through  the  sample 
and  for  about  1  /j. sec  after  shock  exit  from  the  sample.  After  this  time, 
relief  waves  from  the  edges  reach  the  electrodes  and  the  measurement  is  no 
longer  accurate. 

B.  METHOD  OF  MEASUREMENT 

In  the  work  of  Report  I  we  used  a  “ constant- vol tage”  circuit.  For 
the  shots  this  past  year  the  “ constant  -  cur  rent  ”  system  shown  in  Figure  1 
was  used  almost  exclusively. 

In  this  system,  the  power  supply  and  osci 1 loscope  sweep  are  triggered 
5  /J-sec  before  the  shock  reaches  the  specimen.  An  additional  \-fisec  de*ay 
in  the  power  supply  permits  display  of  a  portion  of  the  baseline.  Hence 
the  voltage  drop  across  the  shunt  resistor  is  measured  during  the  early 
portion  of  the  sweep  (the  specimen  resistance  being  effectively  infinite) 
and  the  voltage  drop  across  the  shunt  resistor  and  shocked  specimen  in 
parallel  is  measured  during  the  latter  portion  of  the  sweep.  An  effec¬ 
tively  constant  current  is  obtained  because  the  power  supply  time  constant 
is  many  hundreds  of  microseconds  and  the  ballast  resistance  is  at  least 
fifty  time^.  the  shunt  resistance. 

When  this  system  was  first  used,  it  was  found  necessary  to  use  a 
differential  preamplifier  at  the  oscilloscope  for  shunt  resi stances  < 1 fl; 
that  is,  separate  signal  leads  detected  the  voltage  with  respect  to  ground 
at  each  end  of  the  shunt,  and  the  difference  between  the  signals  was  dis¬ 
played  on  the  oscilloscope.  Even  using  this  system,  we  found  that  the 
signal  level  prior  to  shock  arrival  might  drift  slightly  (usually  decreasing 
in  magnitude)  for  a  shunt  resistance  of  0.02  0,  which  seems  to  indicate 
that  some  contact  resistance  was  varying  with  time.  Since  the  effect  seems 
to  be  unpredictable,  this  resistance  may  depend  on  such  factors  as  the 
quality  of  the  solder  contacts,  although  preliminary  checks  indicate 


*  From  other  work  in  this  laboratory,  the  resistivity  of  this  msterisl  under  shock  compression  is  known 
to  remain  high  at  the  pressures  of  interest. 
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negligible  contact  resistance.  The  maximum  current  used  is  10  amp,  which 
should  cause'  negligible  heating  of  the  resistors  in  the  circuit  for  the 
short  time  of  interest  (^lOyMsec). 

It  was  found  later  that  the  differential  system  was  unnecessary  if  a 
ground  loop  through  the  constant  current  power  supply  was  avoided.  This 
was  accomplished  by  isolating  the  chassis  from  the  common  ground  of  the 
3-’--ire  110-volt  power  and  from  the  delay  unit  by  means  of  a  pulse  trans¬ 
former  at  the  trigger  input. 

The  sample  resistance  R*  is  determined  from  the  shunt  resistance  R_, 
the  signal  prior  to  shock  arrival  ,  and  the  final  signal  V  (measured 
within^O.S  ^ sec  after  the  shock  leaves  the  specimen): 

rx  =  rs[(v0/v)  -  lr1  (i) 

An  example  of  the  type  cf  record  obtained  is  shown  in  Figure  3.  A 
constant-current  system  is  used  with  a  gauge;  a  shunt  resistor  is  not  re¬ 
quired,  of  course,  because  the  initial  resistance  is  finite.  The  Manganin 
resistance  increases  with  pressure  so  that  the  signal  amplitude  increases 
when  the  shock  reaches  the  wire.  A  small  polarization  signal  can  be  seen 

when  the  shock  first  enters  the  insulating  material  (C-7  epoxy)  of  the 
gauge. 

C.  SUBMICROSECOND  MEASUREMENT  OF  LOW  RESISTANCE 

To  gain  a  better  understanding  of  the  behavior  of  the  circuitry, 
several  small-scale  shots  were  fired  with  various  electrode  configurations, 
shunt  resistance  and  inductance,  and  ways  of  handling  electrical  grounds . 

Experiments  in  which  the  shunt  resistance  was  shorted  produced  signals 
which  initially  dropped  very  rapidly,  then  much  more  slowly,  toward  the 
baseline  (see  Figure  4).  Several  microseconds  were  required  before  a  short 
was  indicated.  The  slow  decay  is  consistent  with  themeasured  ~0.  1  in¬ 
ductance  of  the  circuit.  The  initial  drop  is  not  yet  understood,  but 
circuit  inductance  seems  to  be  a  limiting  factor  in  measurements  of  low 
resistance. 

The  time  constant  is  decreased  by  increasing  the  shunt  resistance 
R  ;  for  a  given  R%,  this  means  decreasing  V  (see  Eq.  1).  Since  maximum 
accuracy  is  not  required,  we  can  sacrifice  the  display  of  and  increase 
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(a)  RECORD  FOR  SHOT  NO.  IO,3i3 

UPPER  BEAM -CONDUCTION  SIGNAL 
LOWER  BEAM  -  GAUGE  SIGNAL 


(b)  RECORD  FOR  SHOT  NO.  10,102 
Csl,  Rs-0.l{2 


FIG.  3  TYPICAL  OSCILLOSCOPE  RECORDS  FROM  CONDUCTION 
SHOTS  (to  and  t,  correspond  to  the  shock  entering  and  exiting 
the  sample  respectively.) 


(a)  SHOT  NO  10,090 - 0.02ft  SHUNT 


(b)  SHOT  NO.  10,091 


0.02  ft  SHUNT 


(c)  SHOT  NO.  10,092 - 0.02ft  SHUNT  (d)  SHOT  NO.  10,093 - 0.1ft  SHUNT 

AND  0.02ft  IN  SERIES  WITH  ELECTRODE 

FIG.  4  OSCILLOSCOPE  RECORDS  FROM  SHORTING  SHOTS 


oscilloscope  sensitivity  in  order  to  study  V.  A  Zener  diode  was  used  to 
clip  the  signal  in  order  to  avoid  overdriving  the  oscilloscope,  but  its 
rise  time  was  too  long.  The  diode  was  found  to  be  unnecessary,  however, 
for  reasonable  Rs  values. 

Another  approach  to  the  measurement  of  small  resistances  was  to 
return  to  the  constant- vol tage  circuit  described  in  Report  I.  For  this 
circuit 

Rx  =  RS[(V0/V)  -  1]  (2) 

and  small  values  of  Rx/Rs  mean  V  ~  FQ .  A  Tektronix  Z- type  preamplifier 
was  used  to  expand  the  small  difference  between  V  and  .  The  require¬ 
ment  of  constant  voltage  is  more  stringent  when  the  type  Z  preamplifier 
is  used,  and  the  time  constant  was  consequently  i  r.c  reased  by  i  nc  reas  i  ng  the 
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-Mmnrx  . 


1 

\ 

t 

capacitance.  It  was  necessary  to  place  a  number  of  2-/if  <  apu«  itors  in 
parallel*  to  achieve  the  desired  value,  and  the  inductance  was  thereby 

i 

increased.  No  improvement  over  the  constant-current  circuit  was  achieved. 

* 

D.  PREHEATED  SPECIMENS  {HOT  SHOTS) 

The  design  of  the  hot  shots  (indicated  in  Figure  5)  was  based  on  a 
desire  to  preheat  the  crystal  specimens  to'v'500°C  prior  to  firing  in  a 
configuration  as  close  as  possible  to  that  used  for  room  temperature  shots. 

It  was  found  possible  to  do  this  without  remote  handling  of  the  explosive. 

Both  electrodes  were  heated  externally  by  resistance  wire;  the  tem¬ 
perature  of  each  was  monitored  with  a  separate  thermocouple.  The  bottom 
electrode  (1/4-inch  thick)  was  insulate  from  the  remainder  of  the  driver 
plate  (1/4-inch  thick)  by  a  9-mil  air  gap.  A  Manganin  wire  gauge  rested 
on  a  similar  assembly.  The  aluminum  surfaces  were  polished  to  red  ice  heat 
transfer  and  the  periphery  of  the  driver  plate  was  water-cooled;  nevertheless, 
the  temperature  at  the  center  of  the  driver  plate  exceeded  somewhat  the 
softening  temperature  of  the  explosive  The  explosive  temperature  was  re¬ 
duced  to  an  acceptable  level  by  a  sheet  of  9-mil  Mylar  between  aluminum 
driver  and  explosive.  The  heating  time  to  480°C  was  about  20  minutes. 

Detailed  calculations  of  shocked  states  in  preheated  specimens  are 
presented  in  Appendix  A. 


*  Paper  capacitors  have  been  found  best  for  rspid  release  of  energy;  none  larger  than  2  /if  was  readily 
available . 
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A  SPRING 

B  TOP  HEATER  ASSEMBLY 
C  ASBESTOS 

D  QUARTZ  AND  MICA  INSULATION 
E  HEATER  WINDINGS 
F  SHUNT  RESISTOR 
G  THERMOCOUPLES 


I  CRYSTAL 

J  ALUMINUM  BOTTOM  ELECTRODE 
K  BOTTOM  HEATER  WINDINGS 
L  THERMAL  SHIELD 
M  C-7  EPOXY 

N  RESISTANCE  W'RE  PRESSURE  GAUGE 
0  ALUMINUM 


H  TOP  ELECTRODE 


P  MYLAR 


08  4i00  «« 


FIG.  5  SCHEMATIC  DIAGRAM  OF  HOT  SHOT  ASSEMBLY 
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IV  EXPERIMENTAL  RESULTS 


All  conductivity  shots  fired  during 
in  Table  I.  Except  for  Shots  10,941  and 
made  with  the  constant  -  current  system. 


this  report  period  are  included 
10,942,  all  measurements  were 


The  resistivity  data  from  all  shots  are  plotted  vs.  pressure  in 
Figures  6  and  7,  is.  volume  in  Figures  8  and  9,  and  vs.  3  IT  in  Figures  10, 
11.  and  12.  (The  errors  due  to  the  uncertainties  in  interpreting  the  os¬ 
cilloscope  records  are  discussed  in  Section  V.)  For  the  longitudinal 
geometry,  the  measured  resistance  has  been  converted  to  resistivity  using 
the  compressed  specimen  thickness  and  the  upper  electrode  area. 


The  resistivity  of  NaCl  shocked  to  225  kbar  from  room  temperature 
was  determined  to  be  >  9  *  103  fi-cm.  The  only  shot  for  which  appreciable 
conductivity  was  observed  in  NaCl  was  Shot  No.  10,504,  in  which  the  NaCl 
was  heated  to  480°C  and  then  shocked  to  244  kbar.  The  resulting  resis¬ 
tivity  was  between  186  and  753  fi-cm. 

The  resistivity  of  KI  decreases  to  a  value  of  about  6  Q-cm  at  a 
pressure  of  240  kbar;  a  leveling  off  is  indicated  by  both  the  pressure 
and  volume  plots  (Figures  6  and  8).  Shot  No.  10,939  at  120  kbar  produced 
a  resistivity  of  27  H-cm,  which  is  lower  than  other  results.  The  records 
from  this  shot  do  not  show  any  anomalies,  however,  and  it  has  therefore 
not  been  discarded.  The  initial  temperature  of  all  shots  was  about  25°C. 


For  Cs I  in  the  longitudinal  geometry,  the  resistivity  decreases  from 
M05  fi-cm  at  122  kbar  to  ^6  fi-cm  at  250  kbar.  As  in  the  case  of  KI ,  a 
leveling  off  at  high  pressure  is  indicated  in  Figures  7  and  9. 

It.  should  be  noted  that  the  longitudinal  geometry,  constant -voltage 
data  from  Report  I  do  not  agree  with  the  longitudinal  data  obtained  in 
the  present  work  with  the  constant -current  circuit  for  resistivities 
>  15  fi-cm.  For  this  reason,  two  cons  tant -'  o  1  tage  longitudinal  shots 
(Nos.  10,941  and  10,942)  were  fired  in  the  present  work.  The  results 
agree  well  with  the  present  constant -current  longitudinal  data,  leaving 
a  very  large  discrepancy  between  data  in  this  report  and  those  in  Report  I. 
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Table  1 

SHOT  PARAMETERS  AND  RESISTIVITY  DATA  FOR  ROOi,' 


CRYSTAL 

STATE 

SHOT 

NO. 

Type 

Thickness 
(mm ) 

ELECTRODE 

DIAMETER 

(in.  )° 

SHUNT 

RESISTANCE 

<n> 

LATERAL 

INSUIATION6 

HIGH  EXPLOSIVE C 

un  (  v  f-ii 

MATERIAL  AND 
THICKNESS 

Stress 
(kbar  ) 

Spec 

Vol 

'cm 

9,994 

2.38 

5/8 

25.5 

^"■RJ  ■TUI 

2 "Comp  B-3 

1  '2 

Al ' 

200 

0.3 

10,101  , 

'  ■ 

u 

1.01 

7/8 

0.02 

Bi  fl 

210 

0.3 

10, 206 

■ 

1 

2.18 

7/8 

0.21 

B'f  1 

205 

0.3 

10,312 

1 

2.06 

7/8 

0.20 

Hr  fl  i 

1 

r 

205 

0.3 

9,997 

u 

2.84 

5/8 

7/8 

25.5 

Oil 

2"  9404 

225 

0.3 

10,504 

1 

2.385 

25.5 

Ai  r 

4"  9 

404 

r 

244 

0.3 

10,939 

KI 

1 

1.777 

7/8 

25.5 

Ai  r 

f  i/2"  Comp  B-3 

1 12 

Al 

121 

0.2 

• 

l  +1 "  Bara  to  1 

1.777 

7/8 

25.5 

Ai  r 

2"  TNT 

1/2" 

M 

155 

0.2 

9,998 

1.786 

5/8 

25.5 

Oil 

2 ”  Comp  B-  3 

r 

Al 

190 

0. 1 

■  Wrtr  ■ 

1,013 

7/8 

0.02 

A> 

r 

1/2" 

Al 

205 

0.1 

■  iHvIl 

1.834 

5/8 

0.02 

Air 

206 

0.1 

B  Sjfilijfl 

1.820 

5/8 

0.02 

ni 

207 

0.1 

B  Wf/ */ 

g 

2.03 

7/8 

0.20 

Air 

206 

0.1 

■  J;.’: 

0.988 

7/8 

0.108 

Air 

1 

r 

214 

0. 1 

1.105 

7/8 

0. 106 

Ai  r 

0" 

4. 

9404 

242 

0. 1 

1  BtH 

T 

2.30 

7/8 

3.31 

Air 

2" 

9404 

1 

f 

237 

0.1; 

10,929 

Csl 

1.777 

7/8 

25.5 

Oil 

2"  Baratol 

1/2" 

Al 

122 

0.1 

10,930 

1.777 

7/8 

25.5 

Oil 

• 

T/2"  Cornp  B-3 
+  1*  Baratol 

1/2" 

Al 

154 

0.1 

10,822 

1.676 

7/8 

25.5 

Air 

\ 

T/2"  Comp  B-3 
+  1*  Baratol 

1/2" 

Al 

151 

0.1- 

10,821 

1.83 

7/8 

25.5 

Ai  r 

2" 

TNT 

1" 

Al 

187 

0.1 

10,820 

1.74 

7/8 

25.5 

Air 

2" 

TNT 

1/2" 

Al 

206 

0.1 

9,995 

5/8 

25.5 

Oil 

2"  Comp  B-3 

r 

Al 

230 

0.1 

9,999 

5/8 

0.02 

Oil 

1/2" 

Al 

241 

0.1 

10,097 

7/8 

0.103 

Air 

256 

0.1 

10,306  / 

1.77 

7/8 

0. 17 

251 

0.1 

10,307 

1.51 

7/8J 

0.19 

25  3 1 

0.1 

1.68 

7/8, 

0.23 

2521 

0.;: 

10,308 

1.28 

7/8 

0. 17 

254* 

0.1 

1.29 

7/8, 

0.17 

254* 

0.1. 

10.315 

2.54 

1/4* 

3.40 

247 1 

0. 1 

2.54 

7  '8 

0.20 

247' 

0.1 

10,096 

1.00S 

7/8. 

0.85 

f 

262  . 

0.1 

10,309 

2.34 

1/4* 

3.40 

2" 

9404 

277  * 

0.1 

2.34 

7/8 

0.20 

277  1 

0.1 

10,210 

g 

2.49 

7/8 

0.21 

276  * 

0,1 

2.54 

7/8 

0.18 

276* 

0.1 

10,313  , 

g 

2.49 

7/8 

0.20 

276* 

0.1 

10,314/ 

2.44 

7/8 

3.43 

276  1 

0.1 

10,102 

1.023 

7/8 

0.  i  16 

1 

r 

298 

0.1 

10,100 

1.780 

7/8 

0.110 

3/8 

Al 

(475)/ 

10,503 

h 

2.005 

7/8 

0.195 

2’ 

TNT 

1  '2 

Al 

181 

0.1 

10,931 

h 

1.777 

7/8 

0.265 

2 "  Comp  B-  3 

1 

Al 

230 

0. ! 

10,932 

k 

1.777 

7/8 

0.273 

2"  Comp  B-3 

1 

Al 

0.1 

10,943 

i 

i 

30.0 

1 

11/2"  Comp  B-3 
.  +1*  Baratol 

i 

1/2 

+  1/4 

Al 

Hi  -  D 

118 

0.1 

10,940 

i 

i 

42.8 

2"  TNT 

r  1/2 
+  1/4 

Al 

Hi-1) 

144 

0.1 

10,934 

t 

i 

0.91 

2 "  Comp  B-  3 

1 

r  Mi 
1+1/4 

Al 

Hi-11 

203 

0.1 

10,937 

I 

i 

1.07 

1 

r 

1  « 

4> 

9404 

1 

r  1/2 
L+1/4 

Al 

Hi  -D 

255 

0.1 

10,942 

1.777 

7/8 

25.5 

Air 

2"  Baratol 

1/2 

Al 

129 

0.1 

10,941 

1 

1.777 

7/8 

25.5 

Ai  r 

2 

TNT 

1/2 

AI 

169 

0.1 

Notes : 


0  *%  for  NsCl  and  KI;  A1  for  Cal. 

*  Oil  was  Do*  Corning  200  fluid;  viscosity  =  200  cs. 

C  Initiated  by  a  P-60,  6-inch  plane-wave  generator,  except 
shots  for  which  P-80  planevave  generator  vaa  used. 

^  Average  state  at  tiae  of  measurement;  aee  Table  III. 

*  2024T  A1 . 

'  One 'Ota  If  of  double  shot. 


preheated 

A 


*  Irradiated  90  hours  with  Co6(V  -rays;  2  * 
Crystal  preheated. 

1  Shock  stress  estimated. 

1  See  Sec.  IV  B. 

*  See  Sec .  IV  A. 

*  Transverse  aampfe  1/4  X  1/4  X  0.16  in.,  > 


Table  I 


'ISTIVITY  DATA  FOR  ROOM  TEMPERATURE  AND  HOT  SHOTS 


1 

STATE  IN  SPECIMEN** 

ERIAL  AND 

1 1CKNESS 

Stress 
( kbar  ) 

Sprc  i  f  ic 
Volume  I 

Temperature 

(°K ) 

RESISTIVITY 
( fi  -cm ) 

REMARKS 

(cm  / g ) 

'2* 

Al* 

200 

It”  <  ■ 

1073 

>  5  x  104 

210 

■ft 

1105 

— 

205 

1105 

>  600 

205 

0.3285 

1105 

00  i 

Fired  27  hours  after  2  x  10‘  rad  exposure. 

225 

0.323 

1134 

>  0.9  x  104 

r 

244 

0.326 

2340 

753  -  186 

Crystal  preheated  to  480°C. 

>2- 

Al 

12, 

0.2109 

1645 

28 

'2' 

Al 

155 

0.2015 

2325 

131 

1" 

Al 

0.1940 

2971 

116!$ 

>2' 

Al 

0.1915 

3440 

— 

0.1910 

3550 

207 

0.1908 

3560 

— 

206 

0.1913 

3540 

8-20 

Fired  25  hours  after  2  x  10*  rad  exposure. 

214 

0.1897 

3640 

10  i  5 

242 

0.1848 

4500 

Q+0 

-7 

f 

237 

0.1855 

4540 

2-7 

12' 

Al 

122 

0.1526 

1450 

>  1.6  x  104 

12' 

Al 

154 

0.1484 

2030 

612 

/ 2 

Al 

151 

0.1488 

2075 

160 

r 

Al 

187 

0.1444 

2665 

29 

12' 

Al 

206 

0.1423 

2970 

22 

1" 

Al 

230 

0.1395 

3460 

12  ±  3 

/2" 

Al 

241 

0.1382 

4110 

— 

256 

0.1366 

4260 

5.3^2 

251* 

0.1371 

4240 

1  -  20 

253* 

0.1370 

4240 

2  -  8 

Electrodes  coated  with  silver. 

252* 

0.1371 

4240 

4-14 

254* 

0.1368 

4250 

2-6 

Electrodes  coated  with  graphite. 

254* 

0.1368 

4250 

3  -  20 

247* 

0.1377 

4220 

1.7  -  4.2 

1/4"  electrode  inside  7/8"  OD  guard  ring. 

247* 

0.1377 

4220 

3.3  -  6.7 

262  . 

0.1361 

4375 

5.1  ±  1 

277  1 

0.1344 

5210 

3  -  3.3 

1/4"  electrode  inside  7/8"  OD  guard  ring. 

277  1 

0.1344 

5210 

1.2  -  1.6 

276  1 

0.1345 

5200 

— 

/  Fired  7  hours  after  2  x  10^  rad  exposure.  Trace  appears 

276  1 

0.1345 

5200 

— 

1  normal  for  0,4  psec,  then  goes  rapidly  afcross  base  line. 

276  * 

0.1345 

5200 

2.7  -  4.2 

Fired  30  hours  after  2  x  10?  rad  exposure. 

276  1 

0.1345 

5200 

4 

r 

298 

0.1325 

5600 

6 -8-3 

1/8 

Al 

(475) 1 

— 

— 

(<  0.8) 

Flying  plate  shot. 

/2 

Al 

181 

0.1482 

3620 

1.78 

Crystal  preheated  to  420°C. 

1 

Al 

230 

0.1400 

3780 

<  0.3 
£  0.3 

Crystal  preheated  to  100°C. 

1 

Al 

232 

0.1419 

4540 

Crystal  preheated  to  350°C: 

/2 

/4 

Al 

H;  -D 

118 

0.1530 

1425 

19 

Shot  fired  in  transverse  geometry. 

/2 

/4 

"  Al 

Hi-D 

144 

0.1496 

1820 

(0.2) 

Shot  fired  in  transverse  geometry. 

/2 

/4 

Al 

"  Hi-D 

203 

0.1425 

2965 

0.61  -  1.3 

Shot  fired  in  transverse  geometry. 

1/2 

1/4 

"  Al 
-  Hi-D 

255 

0.1368 

3965 

0. 18 

Shot  fired  in  transverse  geometry. 

1  /2 

■  Al 

129 

0.1515 

1575 

2.3  x  105 

Constant  voltage. 

1/2 

"  Al 

169 

0.1465 

2390 

46 

Constant  voltage. 

W  hours  with  Co^^V-rays ;  2  X  10^  rad  dose, 
heated . 
i  estimated. 

B. 

A. 

sample  1/4  X  1/4  X  0.16  in.;  electrode  width  1/4  in. 


RESISTIVITY - a  -cm 


FIG.  6  RESISTIVITY  vs.  PRESSURE  FOR  KL  POINTS  WITH  ARROWS  INDICATE  LIMITS. 
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RESISTIVITY  vs.  SPECIFIC  VOLUME  FOR  Kl.  POINTS  WITH  ARROWS  INDICATE  LIMITS 


ALDER 

CONSTANT  CURRENT 


FIG.  10  RESISTIVITY  vs.  1/T  FOR  NaCI.  POINTS  WITH  ARROWS  INDICATE 


RESISTIVITY 


Q  -cm 


CONSTANT  VOLTAGE  FROM  REPORT  I 
AL'TSHULER  ET  AL. 


RESISTIVITY  vs.  SPECIFIC  VOLUME  FOR  Csl.  POINTS  WITH  ARROWS  INDICATE  LIMITS 


resistivity 


X 

■< 


ALOER 
HOT  SHOT 

CONSTANT  CURRENT 


RESISTIVITY 


FIG.  (2  RESISTIVITY  vs.  1/T  FOR  Csl.  POINTS  WITH  ARROWS  INDICATE  LIMITS. 
(See  caption  in  Fig.  7) 


25 


The  Cs  I  specimens  used  during  this  report  period  were  of  three  types 
(1)  Harshaw,  nominally  pol ycrysta I  I  ine;  (2)  Harshuw  single  crystals, 
and  (3)  Isotope  single  crystals.  The  specimens  from  which  the  data  ol 
Report  I  were  derived  were  of  type  (1).  The  prese.it  data  for  resistivi¬ 
ties  >  15  il-cm  in  longitudinal  geometry  and  for  all  transverse  shots  were 
obtained  with  type  (2)  crystals  shocked  in  the  <111>  direction.  Thus  the 
discrepancy  may  be  due  to  a  dependence  on  crystallographic  orientation 
and/or  impurity  content. 

In  the  transverse  geometry,  the  resistivities  are  lower  than  for  the 
longitudinal  geometry,  decreasing  to  0.2  ff-cm  at  255  kbar.  Measurements 
in  the  two  geometries  differ  in  two  respects,  viz.,  the  direction  of  cur¬ 
rent  flow  relative  (1)  to  the  crystal  orientation  and  (2)  to  the  direc¬ 
tion  of  shock  propagation.  The  direction  of  current  flow  in  the  transverse 
shots  was  perpendicular  to  the  < 1 1 1 >  direction  but  otherwise  random.  In 
view  of  this  and  of  a  similar  effect  found  in  KI  shocked  in  the  <100> 
direction  (see  Report  I),  a  dependence  of  resistivity  on  direction  of  cur¬ 
rent  flow  relative  to  the  shock  direction  seems  most  likely. 

Several  shots  were  fired  in  which  the  experimental  parameters  were 
varied.  The  results  of  these  shots  and  the  parameters  varied  are  de¬ 
scribed  in  the  following  subsections. 

A.  VARIATION  OF  ELECTRODE  AREA 

For  Shot  No.  10,315,  the  backing  electrode  was  divided  into  a  cen¬ 
tral  0 . 2 50 -i nch -di ame te r  electrode  surrounded  by  a  0 . 87 5  -i nch -0D  x  0.270- 
inch-ID  guard  ring;  the  0.010-inch  annulus  was  filled  with  Teflon.  The 
shunt  resistances  across  the  electrode  and  guard  ring  were  made  inversely 
proportional,  and  the  currents  through  each  directly  proportional,  to 
their  respective  areas.  This  was  to  provide  a  uniform  current  density  in 
the  specimen  (to  the  extent  that  the  resistivity  was  uniform).  The  ad¬ 
vantages  of  this  geometry  are  that  it  eliminates  any  edge  effect  due  to 
stress  relief  at  the  unsupported  edge  of  the  specimen,  and  it  raises  the 
resistance  to  be  measured  by  dt'  reasing  the  effective  specimen  area.  The 
latter  makes  possible  a  large  shunt  resistance  which  should  improve  the 
time  response  of  the  system.  Included  in  this  shot  were  a  second  crystal 
with  a  solid  backing  electrode,  and  a  pi  "ssure  gauge.  ShotNo.  10,309  was 
identical  except  that  it  was  fired  with  9404  rather  than  Comp  B. 


* 

Tins  result  apparently  conflicts  with  the  work  of  Alder  and  Ch r  an — see  Section  V-B. 
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The  data  of  Table  1  show  that  the  electrode  with  a  guard  ring  gave 
essentially  the  same  resistivity  as  a  solid  one,  even  though  the  actual 
resistances  differed  by  an  order  of  magnitude.  The  additional  40  to  50 
kbar  provided  by  9404  compared  with  Comp  B  causes  no  appreciable  change 
in  resistivity.  It  is  worth  noting  that  the  system  using  the  electrode 
plus  guard  ring  should  be  capable  of  indicating  resistivity  a  factor  of 
10  lower  than  that  observed. 

B.  VARIATION  OF  ELECTRODE  MATERIAL 

Shot  Mo.  10,307  was  fired  for  the  purpose  of  comparing  the  resis¬ 
tivity  of  Csl  as  measured  by  using  the  usual  aluminum  electrodes  with 
that  measured  by  using  electrodes  coatea  with  silver  paint.  The  shot 
assembly  consisted  of  two  ervota  Is ,  one  with  regular  electrodes  and  one 
with  coated  electrodes,  plus  a  pressure  gauge.  Shot  No.  10,308  was  iden¬ 
tical,  but  the  coating  was  graphite.  The  results  are  i nc 1 uded  i n  Tab le  I. 
The  coated  and  bare  electrodes  on  a  given  shot  produced  records  which 
were  a  little  different  from  one  another,  but  these  differences  were  no 
greater  than  those  observed  in  other  bare  electric  shots.  The  resistivi¬ 
ties  were  essentially  the  same  for  all  four  crystals. 

C.  VARIATION  OF  SHUNT  RESISTANCE 

To  decrease  the  time  constant  due  to  circuit  inductance,  the  shunt 
resistance  for  Shot  No.  10,314  was  raised  to  3.4  H,  about  10  times  the 
expected  specimen  resistance.  Several  oscilloscopes  were  run  at  rela¬ 
tively  high  sensitivity  so  that  the  signal  after  shock  transit  could  be 
measured,  although  the  initial  signal  due  to  the  shunt  resistor  was  off 
scale.  The  initial  signal  was  observed  on  another  oscilloscope  as  a 
check  on  the  current  and  shunt  resistance.  Again  the  resistivity  was 
several  Q-cm  although  a  value  as  low  as  several  tenths  f)  -  cm  should  have 
been  measurable. 

D.  IRRADIATED  SPECIMENS 

Shots  were  fired  with  crystals  of  NaCl,  K I ,  and  Csl  which  had  been 
irradiated  with  Co60y-rays  2  x  107  rad).  There  appears  to  be  no  sig¬ 
nificant  change  in  the  resistivity  of  NaCl  or  KI ,  but  the  polarization 
signal  in  NaCl  is  modified  somewhat.  On  the  other  hand,  either  irradia¬ 
tion  of  Csl  produced  a  very  large  effect  (of  a  peculiar  nature)  or  the 
shot  was  completely  spurious  for  some  unknown  and  unexpected  reason. 
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K. 


PASS  I  VK  SHOTS  (NO  APPLIED  KIEL!)) 


In  tlio  work  dcsc r i bed  in  He  port  I,  several  passive  shots  were  IT  rod 
w  i  t  h  a  roust  ant  -  vo  I  tagi'  circuit,  and  signals  were  observed  in  every  ruse 
((.si,  kl,  and  NaCI).  TIi i*  nature  of  tin*  signal  varied  greatly  with  speci- 

men  resistivity.  It  sometimes  continued  after  the  shock  exited  from  the 
specimen,  although  the  signal  level  was  small  compared  with  the  conduc¬ 
tivity  signals  and  was  ignored. 

In  the  present  work  it  was  important  to  interpret  correctly  small 
conduction  signals;  therefore,  several  passive  shots  were  fired  (all  with 
Cs  I ) .  Tie  passive  shots  are  listed  in  Table  II  and  a  typical  signal  is 
shown  in  Figure  13(a).  The  record  from  an  active  shot  is  shown  in 
Figure  13(b)  for  comparison.  Not'  that  there  is  a  small  “polarization 
signal”  as  the  shock  traverses  the  specimen  and  that  it  continues  beyond 
the  shock  exit  time  tj — typically  for  1  to  2  ftsec.  The  shots  of  Figure  13 
were  made  with  a  guard  ring  electrode  assembly.  It  is  interesting  that, 
for  both  shets,  the  signal  from  the  center  electrode  exhibits  more  struc¬ 
ture  than  does  that  from  the  guard  ring. 


Table  II 

PASSIVE  SHOT  DATA  FOR  Csl 


SHOT  NO. 

SAMFLE 

THICKNESS 

(  ill  Al  ) 

ELECTRODE 

AREA 

(cm^  ) 

SHUNT 

RESISTANCE 

(ft) 

PEAK 

PRESSURE 

(kbar  ) 

MAXIMUM 

SIGNAL  VOLTAGE'* ' 
(volts) 

DURATION  OF 
SIGNAL  AFTER 
SHOCK  EXIT 
(*isec ) 

10,502 

2.06 

3.88 

0.197 

303 

+0.26,  -0.34 

1.6 

10,817 

1.75 

3.88 

0.27 

275 

+0.30,  -0.55 

1.6 

10,818 

1.55 

3.88 

0.02 

274 

+0.22,  -0.13 

1.6 

10,935(b) 

1.77 

0.316 

3.51 

266 

+0  47,  -0.33 

1.4 

1.77 

3.56 

0.32 

266 

+0.32,  -0.27 

1.1 

10,928 

(c) 

(c) 

1.12 

219 

(d) 

<d) 

a  Positive  voltage  refers  to  voltage  spike  produced  as  ahock  exits  from  specimen.  Nijative 
voltage  refera  to  signal  after  shock  exit  from  specimen  (see  Figure  13a). 

^  Guard  ring  shot;  see  Figure  13a. 

(c  ) 

Transverse  shot;  electrodes  are  aluminum  strips  6.35  mm  wide  and  6.35  mm  apart;  sample 
and  electrodes  are  4  mm  high. 

^  Positive  apike  produced  upon  ahock  entrance  into  sample.  Spike  ia  followed  by  a  gradual 
rise  in  voltage  to  a  maximum  of  +0.17  volt  2)istc  later.  Shock  transit  time  through 
specimen  ia  1.1  usee. 


211 


(a)  RECORD  FOR  SHOT  NO.  10,935 

PASSIVE  SHOT-  NO  APPLIED  FIELD  Csl  AT  266  kbar 
UPPER  BEAM:  CENTER  ELECTRODE,  Rss3.5i2 
LOWER  BEAM:  GUARD  RING  ELECTRODE,  Rs  =  0.32X1 


(b)  RECORD  FOR  SHOT  NO.  10,309 

Csl  GUARD  RING  SHOT  WITH  APPLIED  FIELD  AT  277  kbar 
UPPER  BEAM:  CENTER  ELECTRODE,  Rs=  3.40X1 
LOWER  BEAM:  GUARD  RING  ELECTRODE,  Rs  =  0.302X1 


FIG.  13  COMPARISON  OF  PASSIVE  AND  ACTIVE  SHOTS 
FOR  THE  GUARD  RING  GEOMETRY 
(t1  -  tQ  -  shock  transit  time) 
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One  passive  shot  was  firedwitlithelransver.se  geometry  (Shot  No.  10,938). 
A  pos  i  t  i  ve  -vo  1  tage  spike  was  produced  upon  shock  entrance  into  t!.“  .sample, 
followed  by  a  gradual  rise  in  voltage  to  a  maximum  +0.17  volt  occurring 
2  microseconds  lacer  (shock  transit  time  was  1.1  /xsec) .  The  validity  of 
the  record  is  in  doubt  because  it  was  noisy.  The  signal  did  not  resemble 
the  signals  obtained  with  the  longitudinal  geometry. 

F.  HOT  SHOTS 

The  data  from  four  successful  hot  shots  are  listed  in  Table  1.^ 

The  records  from  Shot  No.  10,503  and  Shot  No.  10,504  are  reproduced  in 
Figures  14a  and  b. 

The  NaCl  record  (Figure  14a)  has  several  strange  features.  In  the 
first  place,  the  characteristic  polarization  signal  produced  as  the  shock 
traverses  the  specimen  is  conspicuously  absent.  The  noise  in  the  record 
beginning  1.2  /x sec  prior  to  the  sharp  rise  in  the  signal  corresponds  to 
the  shock  reaching  the  9 -mi  1  air  gap  separating  the  two  halves  of  the 
driver  plate.  The  mechanism  for  producing  this  signal  is  unknown,  but  we 
have  observed  signals  in  similar  situations  before.  If  the  shock  enters 
the  crystal  at  the  instant  of  the  polarization  signal  produced  in  the 
gauge  ( i.e .,  no  shock  tilt),  the  rapid  increase  in  conductivity  occurs 
about  0.1  /xsec  after  the  shock  enters  the  backing  electrode.  On  other 
shots,  the  onset  of  conduction  corresponds  within  a  few  shakes^  to  the 
shock  reaching  the  electrode.  After  the  initial  rapid  rise,  which  cor¬ 
responds  to  a  resistivity  of  750  O-cm,  the  signal  continues  ’.  rise  less 
rapidly  toward  the  baseline.  The  temperature  of  the  NaCl  in  this  case 
should  cross  the  melting  line  as  the  pressure  is  relieved  to  zero;  there¬ 
fore  this  rapid  drop  in  effective  specimen  resistance  may  be  due  to  melt¬ 
ing  at  the  unsupported  periphery  of  the  crystal.  The  low  resistivity  of 
molten  NaCl,  ~  1  fi-cm  (see  Appendix  B) ,  makes  it  possible  for  a  very 
little  zero  pressure  material  to  effectively  short  out  the  remainder. 

The  Csl  record  is  perhaps  deceptively  good.  There  is  a  positive 
spike  as  the  shock  exits  from  the  specimen.  Fol  lowing  this  spike,  the  volt¬ 
age  across  the  shunt  resistor  remains  nearly  constant  for  about  0.6  /xsec. 

In  fact,  the  voltage  does  ,iot  change  radically  in  the  3.6  /xsec  shown  on  thp 
record.  The  resistivity  measured  after  the  initial  positive  spike  i  s  1 . 8  il-cm . 

* 

The  conduction  signal  is  on  the  order  of  1  volt. 

^Several  shots  were  lost  because  of  electrical  leakage  from  heater  wires  to  electrodes. 

$  1  shake  =  0.01  jxsec. 
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(a)  SHOT  NO.  10,504 - NaCI  SPECIMEN 

PREHEATED  TO  480° C 


(b)  SHOT  NO.  10,503 - Csl  SPECIMEN 

PREHEATED  TO  420°C 

FIG.  14  CONDUCTION  SIGNALS  IN  PREHEATED  SPECIMENS 

(t,  -  t  shock  transit  time) 

i  o 
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In  Shot  No.  10,931,  Cs  I  was  preheated  to  100°C  and  shocked  to 
230  kbar,  and  in  Shot  No.  10,932,  Cs I  was  preheated  to  350°C  and  shocked 
to  232  kbar.  In  both  cases,  the  records  indicate  negligible  voltage 
across  the  shunt  resistor  after  the  shock  exits  from  the  sample,  indicat¬ 
ing  resistivities  of  <  0.3  O-cm. 

G.  TRANSVERSE  GEOMETRY 

Four  Csl  shots  were  fired  in  the  transverse  geometry  with  the  pres¬ 
sures  ranging  from  118  kbar  to  255  kbar.  The  resistivities  determined 
from  these  shots  are  at  least  an  order  of  magnitude  lower  than  the  resis¬ 
tivity  values  from  longitudinal  geometry  shots  (see  Section  IV).  The 
t rans verse -res l s ti vi ty  values  range  from  19  fi-cm  at  118kbar  (Shot  No.  10,943) 
to  0.18  f)-cm  at  255  kbar  (Shot  No.  10,937).  The  value  of  0.2  f)-cm  at 
144  kbar  in  Shot  No.  10,940  does  not  agree  with  the  other  three  trans¬ 
verse  shots.  The  shunt  resistance  used  in  this  shot  was  too  large  so 
that  the  record  indicates  almost  a  short  and  is  hence  less  reliable.  The 
resistivities  quoted  here  and  in  Table  I  are  those  determined  just  as  the 
shock  exits  from  the  specimen. 

The  record  from  a  typical  shot  (No.  10,937)  is  shown  in  Figure  15. 

The  resistance  increases  as  the  shock  traverses  the  specimen  and  then 
levels  off  as  the  shock  exits  and  is  nearly  constant  for  1  psec,  before 
relief  waves  reach  the  electrodes. 


FIG.  15  SHOT  NO.  10,937  Csl  CONDUCTION  SIGNAL 
IN  TRANSVERSE  GEOMETRY.  Rs  -  1.07  Q 

(t,  -  t  =  shock  transit  time) 

I  0 
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For  the  constant  -current  measuring  ci  ji  l,  the  sample  res  i  s  la  n<  <•  H 
was  given  b\  Eq .  1  as 

=  RjU'  /V)  ~  l]'1  (I) 


The  sample  resistance  can  also  be  written  in  terms  of  the  sample  resis¬ 
tivity  p,  the  cross -sections  1  area  of  the  sample  A,  and  the  length  of  the 
current  path  L,  as 


Rx  -  pL/A  .  (3) 

If  these  two  expressions  are  equated  and  ratio  V /V  is  solved  for,  the  re¬ 
sult  is 

V/V0  =  (1  +  RA/pL)-'1  .  (4) 

For  the  transverse  geometry,  the  cross-sectional  area  of  the  conduct¬ 
ing  region  varies  with  time  as  the  shock  wave  traverses  the  specimen,  and 
is  given  by 


A  *  w(U  -  u)t 


(5) 


where  w  is  the  width  of  the  specimen,  11  is  the  shock  velocity, 
particle  velocity  behind  the  shock,  and  t  is  the  time  In  all 
transverse-geometry  shots,  w  =  L.  Substituting  Eq .  5  into  Eq 


u  is  the 
the 

4  yields 


V/V  0  =  [1  +  Rs(U  ~  u)t/pVX 


(6) 


Figure  16  shows  a  plot  of  Eq .  6  for  the  parameters  of  Shot  No.  10,937 
compared  with  the  observed  value  of  F/V0  taken  from  the  record  shown  in 
Figure  15.  The  value  of  o  used  for  Figure  16  is  that  determined  from  the 
record  of  Figure  15  at  the  time  the  shock  leaves  the  crystal. 

The  two  curves  in  Figure  16  do  not  have  exactly  the  same  shape. 

This  result  is  possibly  due  to  attenuation  of  the  shock  in  the  specimen. 
An  alternative  way  of  determining  the  resistivity  would  be  to  match  the 
initial  slopes  of  the  two  curves.  If  this  is  done,  the  resulting  value 
of  the  resistivity  is  greater  by  a  factor  of  two  than  the  resistivity 
determined  at  the  time  of  shock  exit. 
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FIG.  16  PLOT  OF  CONDUCTION  SIGNAL  vs.  TIME  FOR  THE  CASE  OF  TRANSVERSE 
GEOMETRY.  CALCULATED  AND  OBSERVED  VALUES  ARE  SHOWN. 
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V  DISCUSSION  OF  DATA 


A.  PRECISION 

In  Figure  17,  oscilloscope  records  are  shown  for  shots  at  high  pres¬ 
sure  with  various  values  of  shunt  resistance.  It  seems  clear  that  a 
large  increase  in  electrical  conductivity  is  indeed  produced  by  a  strong 
shock.  The  highest  conductivity  measured  in  the  longitudinal  configura¬ 
tion  is  near  the  upper  limit  of  measurement.  However,  as  discussed  in 
Section  IV,  an  order  of  magnitude  higher  conductivity  should  have  been 
measurable  using  a  high  shunt  resistance  or  a  guard  ring  type  backing 
electrode.  The  shape  of  the  signals,  however,  is  poorly  understood,  and 
therefore  the  interpretation  of  oscilloscope  records  is  quite  subjective. 
The  signal  amplitude  may  vary  considerably  during  the  first  microsecond, 
particularly  if  the  shunt  resistance  is  small  and  approximately  equal  to 
the  specimen  resistance.  Some  light  is  shed  on  this  by  the  passive  shots 
discussed  above,  which  indicate  that  at  least  part  of  the  observed  struc¬ 
ture  in  the  signal  during  the  first  /x sec  following  shock  transit  is 
generated  by  the  specimen. 

The  question  naturally  arises  as  to  what  influence  this  generated 
voltage  has  on  the  measurement  of  the  specimen  resistivity.  In  an 
attempt  to  answer  this  question,  the  results  of  active  and  passive 
shots  (see  Figure  13)  for  similar  explosive  systems  and  shunt  resistors 
were  superimposed  on  one  another.  This  comparison  suggests  that  the 
resistivity  values  which  have  been  previously  reported  are  not  grossly 
in  error. 

There  are  several  things  which  make  this  approach  uncertain,  however 
The  crystal  thickness  of  the  active  shot  was  not  identical  to  that  of  the 
passive  shot  used  for  comparison,  nor  were  the  shunt  resistors  of  the  com 
pared  shots  identical.  In  a  majority  of  the  shots  compared,  the  net 
voltage  became  positive*  within  1  ,u sec  after  the  shock  exited  from  the 
specimen.  Possibly  more  important  is  the  fact  that  the  conductivity 


* 

A  conduction  signal  is  negative. 
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(a)  SHOT  NO  10,097 


0.1  ft  SHUNT 


(b)  SHOT  NO.  10,096 - 085ft  SHUNT 


(c)  SHOT  NO.  10,314 - 343ft  SHUNT 


FIG.  17  COMPARISON  OF  RECORDS  AT  HIGH  PRESSURE  FOR  DIFFERENT  Rs 
(t  -  t  -  shock  transit  time) 
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record  is  obtained  with  an  electric  field  applied  across  the  sample, 
whereas  the  record  of  the  voltage  generated  by  the  specimen  is  obtained 
with  no  applied  field.  It  is  not  known  how  the  presence  of  an  electric 
field  will  influence  the  voltage  generated  by  the  crystal.*  It  is  par¬ 
ticularly  noteworthy  in  the  passive  records  that  a  sharp  positive  spike 
is  observed  as  the  shock  reaches  the  backing  electrode.  This  justifies 
our  neglect  of  this  spike  in  interpreting  earlier  conductivity  records. 

The  thermodynamic  state  of  the  specimen  is  also  uncertain  because 
of  attenuation  of  tie  shock  wave  and  the  small  impedance  mi smatch * wi th 
the  backing  electrode.  The  temperature  has  a  further  uncertainty  because 
it  must  be  calculated  using  approximate  thermodynamic  parameters.  At  the 
beginning  of  this  study,  precise  knowledge  of  the  pressure  was  considered 
of  secondary  importance,  and  pressures  reported  in  Report  I  were  those 
measured  with  similar  explosive  assemblies  on  other  projects.  At  the 
same  time  that  the  extreme  pressure  dependence  of  conductivity  in  Csl 
was  being  discovered,  a  resistance  wire  pressure  gauge  was  being  developed 
by  another  group  in  the  Poulter  Research  Laboratories.6  This  gauge,  used 
routinely  on  shots  during  the  past  year,  not  only  indicates  peak  pressure 
but  is  believed  to  follow  the  subsequent  pressure  decay — how  faithfully 
and  for  how  long  are  subjects  of  current  investigations.  The  duration 
of  fidelity  is  primarily  determined  by  gauge  dimensions  and  the  gauge 
used  in  this  work  should  indicate  shock  decay  at  least  semi-quantitative ly 
during  the  first  microsecond.  We  have  used  these  data  to  correct  indicated 
peak  pressures  to  approximate  average  pressures  in  the  specimen.  This 
correction  is  indicated  in  Table  III  for  the  various  explosive  assemblies 
used.  It  should  be  noted  that  the  data  from  Report  I  which  are  used  in 
this  report  have  also  been  corrected. 

The  pressure  correction,  although  approximate,  is  obtained  directly 
from  gauge  records.  The  corresponding  corrections  to  the  volume  and 
temperature  are  not  so  straightforward.  They  have  been  estimated  by 
assuming  that  the  pressure  behind  the  shock  releases  along  the  Hugoniot, 
i.e.,  that  the  isentrope  is  approximated  by  the  Hugoniot  for  the  pur¬ 
pose  of  adjusting  the  volume  and  temperature.  The  result  is  that  the 
volume  and  temperature  corrections  are  underestimated,  as  shown  by  an 
example  in  Appendix  C. 


In  terms  of  the  field  required  to  produce  a  volume  polarization  equivalent  to  that  produced  by  the  shock, 
the  field  applied  in  these  experiments  should  produce  a  negligible  effect. 
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Table  111 


CORRECTIONS  FOR  SHOCK  DECAY 


EXFLOSIVE 

ASSEMBLY 

PEAK  STRESS 

IN  DRIVER* ** 

(kbar ) 

SPECIMEN 

PEAK  STRESS 

IN  SPECIMEN 

(kbar) 

AVERAGE  STRESS 
IN  SPECIMEN*1 
(kbar) 

P- 60  +2"  Baratol 
♦  1/2"  A1 

150 

Csl 

127 

125 

P- 60  +  1/2“  Comp  R-3 

183 

kl 

129 

121 

♦  1"  Baratol 

♦  1/2"  Ai 

Csl 

161 

154 

P-60  +  2”  TNT 
v  1"  Al 

221 

Csl 

198 

187 

P-60  •»  2"  TNT 

223 

KI 

169 

155 

♦  1/2"  Al 

Csl 

214 

205 

P-60  +  2 “  Comp  B-3 

267 

KI 

204 

190 

+  1"  Al 

Csl 

244 

230 

P- 60  +2"  Comp  B-3 

305 

NaCl 

231 

210 

+  1/2"  Al 

KI 

235 

206 

Csl 

283 

250 

P-60  +  2"  9404 

352 

NaCl 

261 

226 

+  1/2"  Al 

KI 

275 

238 

Csl 

330 

277 

P-80  +  2"  TVT 
+  1/2"  Al 

221 

CsI(Hot) 

18« 

181 

P-80  +  4"  9404 
+  1/2"  Al 

382 

NaCl (Hot) 

279 

244 

P-60  +  2'  9404 
+  1/2"  Al 
+  1/4”  Hi-D 

314(in  Hi-D) 

Csl 

267 

255 

P-60  +2"  Comp  B-3 
+  1/2"  Al 
+  1/4"  Hi -D 

257 ( in  Hi-D) 

Csl 

214 

203 

P-60  +  2”  TNT 
+  1/2"  Al 
+  1/4"  Hi-D 

1 84( in  Hi-D) 

Csl 

146 

144 

P-60  +  1/2"  Comp  B-3 
+  1"  Baratol 
+  1/2”  Al 
+  1/4"  Hi-D 

1 56 (in  Hi-D) 

Csl 

121 

118 

*  Except  for  hot  shots,  the  pesk  stress  in  driver  is  the  average  of  several  shots, 

**  At  time  of  measurement.  This  is  tsken  to  be  0.5  *isec  after  shock  reaches  backing  electrode.  The 
pressure  decresse  depends  also  on  specimen  thickness,  the  listed  values  being  representative  ones. 


38 


Hie  Hugoiiiols  for  the  specimen  materials  are  from  the  work  of 
Christian8  and  are  ->hown  in  figure  18,  along  with  curves  for  aluminum 
and  magnesium,9  the  electrode  materials.  The  tabulated  specimen  pressures 
are  for  the  initial  transit  of  the  shock  through  the  specimen;  no  account 
is  taken  of  the  small  mismatch  at  the  specimen-backing  electrode  interface. 

Calculated  shock  temperatures  according  to  Christian8  and  Al’tshuler 
et  a l .  10  (NaCl  only)  are  plotted  in  Figure  19.  Christian  used  the  simple 
technique  given  by  Walsh  and  Christian11  in  which  Griineisen’s  parameter  / 
is  assumed  to  be  proportional  to  the  volume  V.  The  more  elegant  approach 
of  Al’tshuler  et  al.  12  predicts  y  decreasing  less  rapidly  than  V  for  Csl. 
Their  calculated  temperature  data  at  low  pressures  look  abnormally  high, 
however,  and  the  temperatures  used  in  subsequent  discussion  will  be  those 
according  to  Christian.* 

B.  COMPARISON  WITH  OTHER  WORKS 

The  work  of  Alder  and  Christian1  which  prompted  this  study  was  pre¬ 
sented  more  quantitatively  in  a  later  paper  by  Alder. 13  These  later  data 
for  NaCl,  KI,  and  Csl  are  included  in  Figures  6-12;  the  agreement  with 
our  data  is  only  fair. 

In  the  case  of  Csl,  Alder’s  260-kbar  point  agrees  with  our  transverse 
data,  and  his  190-kbar  point  agrees  reasonably  well  with  our  longitudinal 
data.  Alder  and  Christian  reported  that  they  observed  no  difference  between 
longitudinal  and  transverse  geometries,  although  the  details  of  their 
experiments  are  sketchy.  For  KI,  the  agreement  is  fair;  Alder’s  point 
falls  within  the  scatter  of  our  data. 

The  NaCl  data  disagree  with  both  Alder  and  Al’tshuler  et  al.;  our 
resistivities  are  at  least  an  order  of  magnitude  higher  than  theirs. 
Al’tshuler  e  t  a l .  indicate  that  the  scatter  in  their  data  results  from 
using  crystals  from  different  sources,  and  hence  of  presumably  different 
puri ties. 


We  have  checked  these  temperatures  and  obtained  reasonable  agreement  with  Christian's. 
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FIG.  18  PRESSURE  vs.  PARTICLE  VELOCITY  FOR  NoCI,  KI,  Csl,  Mg,  AND  Al 
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FIG.  19  PRESSURE  vs.  SHOCK  TEMPERATURE  FOR  NaCl,  KI,  AND  Csl 


VI  THEORETICAL  CONSIDERATIONS 


We  have  not  yet  determined  what  the  conduction  mechanism  is;  in  the 
following  discussion  we  consider  the  experimental  results  within  the 
framework  of  two  different  hypotheses. 

A.  ELECTRONIC  CONDUCTION 

Alder 13  he  1 ie ves  the  current  carriers  are  electrons  and  originally 
labeled  the  effect  a  “metallic  transition.”  He  suggested  that  the  band 
gap  is  being  decreased  bv  pressure  to  such  an  extent  that  electrons  are 
thermally  excited  into  the  conduction  band.  Considering  the  compressed 
specimen  to  be  an  intrinsic  semiconductor  and  making  assumptions  about 
the  mobility  to  arrive  at  a  pre -exponential  factor,  he  derived13 


,0(0  -  cm)  =  ( 5 . 45/T)eCg/2  * T  (7) 

where  Eg  is  the  band  gap.  The  values  of  Eg  given  by  Alder,  corresponding 
to  measured  p  and  calculated  T,  a^e  shown  in  Table  IV  with  similar  values 
calculated  from  our  data.  The  values  of  Eg  at  1  atmosphere  (Eg°;  are  also 
given  and,  because  they  are  determined  optically,  should  be  upper  limits 
for  the  thermal  band  gaps. 

The  single  point  of  Alder  md  Christian  for  NaCl  agrees  with  our 
calculated  band  gaps.  The  data  indicate  that  in  the  200-  to  244-kbar 
range  the  calculated  tnermal  band  gap  is  3  to  5  ev  compared  to  an  optical 
band  gap  of  8.4  ev  at  atmospheric  pressure  and  ambient  temperature. 

Alder  and  Christian’s  value  of  4.1  ev  at  160  kbar  for  KI  agrees  very 
well  with  our  value  of  4.4  ev  at  155  kbar.  KI  undergoes  a  phase  change 
at  a  pressure  much  less  than  160  kbar  ar.d  it  is  not  known  what  band  gap 
to  associate  with  the  high  pressure  phase.  However,  our  data  combined 
with  Eq.  7  indicate  an  increase  of  about  2.2  ev  in  the  band  gap  -s  KI  is 
further  compressed  (and  the  temperature  increased)  to  242  kbar. 
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BAND  GAP  CALCULATIONS 
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for  (...-  I  ,  t  lie  band  gaps  calculated  I  i  oia  the  longitudinal  shot  data 
•lie  roiisi  strut  ly  (rrcatc  r  tli.m  those  calculated  i  rom  tin-  transverse  shot. 

;i  t  ii .  lhe  2b0-kbar  point  ol  Aider  and  Christian  agrees  with  oil  r  trans- 
xerse  data,  anil  l  lie  i  r  1 QQ  kbar  point  agrees  with  our  longitudinal  data. 

As  previously  discussed,  Alder  and  Christian  report  no  difference  in  their 
data  for  longitudinal  and  transverse  geometries. 

Preheating  the  specimen  was  a  first  attempt  to  separate  the  volume 
and  temperature  dependence  of  Lhe  conductivity.  The  Cs I  shot  with  TNT  was 
designed  to  reach  approx  mate ly  the  same  final  temperature  as  a  room  tem¬ 
perature  Comp  B  shot,  but  at  a  lower  pressure  (and  larger  volume).  (The 
method  of  calculating  the  Hugoniot  ol  the  preheated  material  is  the  subject 
of  Appendix  A. )  The  calculated  states  and  measured  resistivities  are 
shown  in  Table  I.  Assuming  Eq .  7  to  be  valid.  Eg  is  given  in  terms  of 
the  resistivities  at  two  different  temperatures  but  the  same  volume  by 


Eg  = 


2k  in  (, Pjl/P2T2 ) 


(8) 


Using  Eq.  8  and  the  values  in  Table  I  gives  Eg  =  4.2  ev  for  Csl.  Using 
Eq.  7  separately  for  each  point  yields  values  of  Eg  of  4. 3-4. 4  ev,  but 
the  use  of  Eq.  8  seems  most  valid  because  the  questionable  magnitude  oi 
the  p re- exponent i a  1  term  is  not  involved. 

Following  the  initial  apparent  decrease  in  the  band  gap,  the  data 
from  both  geometries  indicate,  as  for  KI,  an  increase  in  band  gap  as  the 
shock  strength  is  increased.  Three  obvious  causes  of  such  a  trend  are: 

(1)  the  measured  values  of  resistivity  are  too  high  at  high  pressures; 

( _ )  The  calculated  temperatures  are  too  high;  and  (3)  Eq .  7  is  not  of  the 
proper  form.  We  have  convinced  ourselves  that  our  experimental  technique 
is  capable  of  measuring  resistivities  an  order  of  magnitude  lower  than 
those  observed.  Although  there  is  considerable  uncertainty  in  the  calcula  - 
tion  of  the  shock  temperature,  there  is  no  reason  to  believe  that  the 
calcula ted  values  used  here  are  excessive;  indeed,  the  temperatures  pub¬ 
lished  by  Al'tshuler  et  al.  are  considerably  higher.  Equation  7  certainly 
is  vulnerable  to  criticism.  The  principal  assumptions  involved  in  its 
derivation  are  that  the  specimen  behaves  as  an  intrinsic  semiconductor, 
and  that  the  mobility  varies  as  T  1  ^  2 .  In  Appendix  D  an  alternative 
expression  is  derived,  which  retains  the  former  assumption  but  replaces 
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the  latte  I 
r  1  iikmi  t 
dclonill  lie 


iiv  a  T~ '■  dcpcndeii  e  e .  I  li  i  >  more  closely  approximates  tli 
;|  I  <lata  lor  polar  semiconductors,  surli  data  being  used  a  1 
a  pre- exponent  i  a  I  constant  (see  Appendix  I)).  The  result 


r 


so 
i  s 


to 


r  =  (1.2  x  10'5)7'e£*/2fcT 


(9) 


Band  gaps  calculated  from  Eq.  9  are  also  included  in  Table  IV.  The  band 
gaps  are  lower  than  those  calculated  from  Eq.  7,  and  the  change  with  in¬ 
creasing  pressure  is  not  as  great. 

Of  course  no  possible  effects  of  shock  compression  on  the  crystal 
lattice — e.g.,  the  influence  on  the  mobility  of  shock-produced  point 
defects,  dislocations,  perturbed  phonon  spectrum — are  included  in  Eqs.  7 
and  9.  Again,  if  impurity  conduction  is  important,  Eqs.  7  and  9  are  not 
appl icable. 

The  data  from  the  longitudinal  and  transverse  geometries  lie  on 
separate  curves,  the  transverse  geometry  indicating  the  lower  resistivi¬ 
ties.  Two  possible  mechanisms  are  suggested  for  a  directional  dependence. 
To  the  extent  of  its  ability  to  support  a  shear  stress,  the  cubic  symmetry 
of  the  crystal  is  perturbed  under  compression  by  a  plane  shock  wave;  hence 
the  conductivity  need  not  remain  a  scalar  quantity.  Also,  the  mobility  is 
influenced  by  lattice  imperfections;  the  distribution  of  dislocations,  for 
example,  is  not  expected  to  be  isotropic. 

This  interpretation  in  terms  of  decrease  in  band  gap  with  latti’e 
compression  is  in  conflict  with  the  theoretical  work  of  Flower  end  March.5 
They  calculated  the  band  gap  of  Cs I  as  a  function  of  isotropic  compression 
of  the  lattice  to  250  kbar  and  concluded  that  any  change  was  negligible. 

B.  IONIC  CONDUCTION 

X) 

Al’tshuler  et  al.  suggested  that  the  conductivity  they  observed  in 
NaCl  at  very  high  pressures  was  ionic  in  nature.  This  conclusion  was 
based  on  a  linear  inp  vs.  1 /T  plot  giving  an  activation  energy  of  1.2  ev, 

a  magnitude  similar  to  that  observed  for  ionic  conductivity  at  1  atmo- 

•  1  *  ... 

sphere.  *■*  They  gave  no  discussion  in  terms  of  a  pressure  or  volume  dependence. 

Our  data  for  NaCl  are  too  few  and  at  too  low  a  pressure  to  confirm 
Al’tshuler’s  suggestion.  Figure  10  shows  a  plot  of  Al’tshuler’s  data 
(the  dashed  lines  are  his).  The  resistivity  determined  for  our  hot  shot 
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is  an  order  of  magnitude  greater  Ilian  his  at  the  same  temperature.  it  is 
tempting  to  conclude  that  this  fact  indicates  a  volume  dependence  of  the 
resistivity  —  the  volume  for  our  fiot  shot  being  greater  than  the  volume 
lor  A1 ’ tshuler ’ s  room  temperature  shot.  This  is  unwarranted,  however, 
because  our  lower  pressure  room  temperature  shots  also  exhibit  higher 
resistivities  than  those  of  the  Russians. 

for  both  KI  and  Cs ! ,  our  data  on  a  log  (j  vs.  1/7  plot  can  he  roughly 
fitted,  at  the  higher  temperatures  (and  higher  pressures),  by  a  straight 
line  (see  Figures  11  and  12).  The  activation  energy  determined  from  the 
slope  is  1.6  ev  for  KI  and  1.1  for  Csl.  The  data  lor  Cs I  are  compared 
with  zero  pressure  data  in  Figure  20. 

At  1  atmosphere,  the  conductivity  of  the  alkali  halides  is  by  dif¬ 
fusion  of  ions  (usually  positive  at  lower  temperatures,  with  the  contribu¬ 
tion  of  negative  ions  increasing  with  temperature)  via  a  vacancy  mechanism. 
The  magnitude  of  conductivity  is  determined  by  the  concentration  of  vacan¬ 
cies  and  by  their  mobility.  If  the  crystal  is  maintained  in  the  thermo¬ 
dynamic  equilibrium  the  conductivity  will  be  intrinsic  at  high  temperature, 
i.e.,  independent  of  the  impurity  content.  Below  some  fairly  well  defined 
temperature  (which  increases  with  impurity  content),  the  concentration  of 
vacancies  is  determined  primarily  by  the  presence  of  divalent  impurities; 
this  region  is  called  extrinsic. 

Although  the  equilibrium  temperature  dependence  of  conductivity  for 
the  materials  used  in  this  study  are  relatively  well  known  (see  Figure  20), 
the  pressure  dependence  is  not.  The  only  measurements  on  ionic  solids 
above  1  kbar  are  those  of  Pierce16  to  9  kbar  on  the  chlorides  of  Na ,  K, 
and  Rb.  These  measurements  are  in  the  extrinsic  range;  furthermore,  the 
pressure  dependence  is  characterized  by  an  activation  volume  (see  Appendix  B 
for  discussion  of  ionic  conductivity)  which  shows  no  systematic  pressure 
dependence  in  the  low-pressure  range  available  to  Pierce,  but  which  must 
certainly  be  pressure  -  dependent .  Hence  Pierce's  data  merely  provide  a 
starting  point  for  estimating  the  pressure  dependence  at  the  pressures 
of  interest.  Such  an  estimate  is  discussed  in  Appendix  B.  The  result 
is  that  one  cannot  say  much  more  than  that  for  weak  shocks  (low  tempera¬ 
tures),  the  pressure  dependence  overshadows  the  temperature  dependence, 
and  that  the  opposite  is  true  for  strong  shocks  (high  temperatures). 

Thus  the  conductivity  should  first  decrease  and  then  increase  as  the 
shock  strength  is  increased.  One  might  expect,  then,  that  the  usual 
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FIG.  20  PLOT  OF  Inp/T  vs.  1/T  FOR  Csl  SHOWING  HIGH  PRESSURE 
DATA  IN  RELATION  TO  ZERO  PRESSURE  DATA.15 
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lnrT  r.s.  1  T  plot  (.’  =  conductivity)  lor  shock  experiments  would  have  a 
larger  slop*1  at  !i  1  gh  temperatures  (corresponding  to  strong  shocks;  than 
a  similar  plot  at  P  =  0,  and  hence  would  have  an  abnormally  high  apparent 
value  of  activation  energy. 

Another  approach  to  the  problem  of  pressure  dependence  is  to  assume 
that,  at  equilibrium,  the  conduc t i v i ty  i  s de t e rmi ned  by  T/Tm  ,  where  Tm 
is  the  melting  point  —  i.e.,  a  corresponding  states  approach.  Again, 
because  the  ratio  of  Hugoniot  temperature  to  melting  temperature  first 
decreases  and  then  increases  as  shock  strength  is  increased,  the  con¬ 
ductivity  would  do  likewise. 

There  is  still  another  argument.  The  pressure  dependence  of  the 
conductivity  of  molten  salts  is  probably  small  (see  Appendix  B).  It 
seems  reasonable  that  the  conductivity  change  on  melting  (increases  by 
MO3  at  P  =  0)17  should  be  relatively  insensitive  to  pressure.  Hence  the 
pressure  dependence  of  conductivity  in  the  solid,  near  the  melting  point, 
should  be  small. 

The  reason  for  this  emphasis  on  the  region  near  the  melting  curve  is 
that  the  temperatures  behind  the  strong  shock  waves  producing  high  conduc¬ 
tivity  in  KI  and  Cs 1  are  probably  near  the  melting  point. 

The  only  pertinent  data  on  melting  curves  of  alkali  halides  are  those 
of  Clark,  who  made  measurements  to  25  kbar  on  the  sodium  halides  and  the 
alkali  chlorides.18  There  is  little  choice  but  to  use  the  two-parameter 
Simon  equation,  with  Clark’s  values  of  the  constants,*  to  extrapolate  an 
order  of  magnitude  in  the  pressure.  The  results  of  such  extrapolations 
are  shown  in  Figure  21  along  with  the  calculated  shock  temperatures.  The 
uncertainties  in  the  melting  curve  for  CsCl-II  are  those  indicated  by 
Clark;  the  uncertainty  is  similar  for  NaCl  and  Nal  ana  larger  for  KCl-II 
(one  parameter  was  assumed  because  phase  II  was  present  only  over  a  5*kbar 
region).  It  must  be  concluded  that ,  if  there  is  time  for  melting  to  occur, 
there  is  a  good  chance  that  our  specimens  are  in  the  molten  state  for  the 
strongest  shocks.  Available  evidence  indicates  that  melting  can  occur, 
at  least  during  pressure  release,  in  the  brief  duration  of  a  shock  experi¬ 
ment.  This  evidence — admittedly  indirect  —  is  the  failure  of  measured 


Recorrecting  Clark'?  thermocouple  measurements  of  temperature  for  the  effect  of  pressure,  using  the 
recent  work  of  Bundy  and  Strong,  *has  negligible  effect  on  the  melting  curves  in  comparison  with  the 
uncertainties  in  determining  the  parameters  in  the  Simon  equation# 
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FIG.  21  HUGONIOT  TEMPERATURES  (solid  lines)  ACCORDING  TO 
CHRISTIAN8  (C)  AND  AL'TSHULER  etal,10'12  (A). 
Melting  curves  (dashed  lines)  extrapolated  from  25  kbar  data 
af  Clark. 1  8  The  suffix  H  refers  ta  the  polymorph  having  the 
CsCI  structure.  Presumably  the  iodides  and  chlorides  have 
nearly  identical  melting  curves  (true  for  Na). 


f  it'f  -  sur  1  ace  velocity  to  equal  twice  the  particle  velocity  (as  determined 
by  the  impedance  mismatch  technique)  for  low  melting-point  metals9  and 
for  the  alkali  halides8  themselves. 

Yaffe  and  Van  Ar tsda 1 en 20  ha \ e  measured  the  conductivities  of  several 
molten  salts  at  1  atmosphere.  They  found  values  of  the  order  of  1  (O-cmT1 
with  a  very  small  temperature  dependence,  which  would  result  in  activation 
energies  of  0. 1-0.2  ev  for  NaCl ,  Cs I ,  and  KI,  one  order  of  magnitude  lower 
than  for  the  solid  state.* 

In  summary,  we  might  expect  the  following  behavior  of  ionic  conductivity- 
assuming  thermodynamic  equilibrium.  The  temperature  increases  slowly  for 
weak  shocks;  we  expect  the  pressure  effect  to  be  dominant  and  the  conduc¬ 
tivity  to  decrease  with  increasing  shock  strength.  At  higher  shock 
strengths  the  temperature  begins  to  increase  so  rapidly  that  the  pressure 
effect  becomes  minor  and  the  downward  trend  of  conductivity  is  reversed. 

In  particular,  near  the  melting  line  we  expect  conductivities  of  the  order 
of  10“3  to  10~4  (0  ■cm)’1 ,  as  at  1  atmosphere;  likewise,  if  melting  occurs 
we  expect  ~1  (fi'em)’1  and  increasing  the  shock  strength  should  have  little 
e  f  feet . 

The  heated  Cs I  shot  is  consistent  with  the  melting  hypothesis.  It 
may  be  that  a  specimen  initially  at  room  temperature  melts  when  shocked 
to  270  kbar.  If  this  is  so,  it  will  certainly  melt  when  shocked  to 
approximately  the  same  temperature  at  a  lower  pressure  ('N'180  kbar).  The 
corresponding  resistivities  are  5  and  1.8^-cm,  respectively.  The  data 
can  also  be  used  to  compare  resistivities  at  different  temperatures  but 
the  same  volume.  Using  Eq.  B-3,  the  activation  energy  is  given  by 


k  In  (PjT2/  p2T1 ) 


(10) 


the  ionic  analog  of  Eq.  8.  From  the  data  of  Table  I  and  Eq.  10, 
i\G  =  2.5  ev  (at  a  volume  corresponding  to  154  kbar  on  the  . oom  tempera¬ 
ture  Hugoni ot )  which  is  double  the  value  of  1.25  ev  .ound  for  Cs I  in  the 
intrinsic  region  at  zero  pressure.  However,  if  melting  is  occurring, 


* 


It  haa  recently  been  obaer»ed  that  the  co 
temperature  well  above  the  melting  point. 


^ductivity  of  molten  salta  may  decrease  with  increaaing 
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most  o  I  l  In*  drop  in  resistivity  run  In*  nsi  lilted  to  it.  Add  to  this  a 
possible  pressure- produced  increase  in  resistivity  contributing  to  the 
low  t  en.pe  ra  t  ure  point,  and  one  can  easily  rationalize  the  larger  value 
o  1  What  is  not  so  easy  to  rationalize,  if  melting  is  occurring,  is 

the  discrepancy  betwi  *n  the  longitudinal  and  transverse  experiments  and 
the  discrepancy  between  experiments  on  different  specimens,  which  we  have 
tried  to  ascribe  to  variations  in  impurity  content  or  a  dependence  on 
crystallographic  orientation. 

Implied  in  this  discussion  is  the  as  umption  that  the  pressure  ’  j 
hydrostatic  and  the  process  of  compression  is  quas  i -s ta t  i  c .  Let  us  con¬ 
sider  the  consequences  of  admitting  that  we  have  shocked  the  crystal.  In 
the  first  place,  dislocations  are  presumably  being  produced  and  driven  at 
high  velocity  and  they  may  be  leaving  vacancies  and  interstitials  behind. 
Both  the  dislocations  and  the  point  imperfections  can  contribute  to  the 
conductivity.  In  the  second  place,  our  measurements  are  made  during  the 
first  /xsec  behind  the  shock  front,  which  seems  a  rather  short  time  in 
which  to  consider  a  process  quas i -s ta t i c .  However,  the  jump  frequencies 
become  very  high  at  high  temperatures,  and  the  dislocations  provide  a 
large  concentration  of  sources  and  sinks  for  point  imperfections,  so 
that  thermodynamic  equilibrium  may  be  possible  in  the  shocked  solid  near 
the  melting  line.  This  will  further  enhance  the  increase  of  conductivity 
to  be  expected  as  the  shock  strength  (and  hence  temperature  and,  specifi¬ 
cally,  T/T  )  is  increased. 

C.  OTHER  POSSIBLE  MECHANISMS 

Other  conduction  mechanisms  include  the  following  poss  bilities: 

1.  Both  electrons  and  ions  may  be  contributing  to  the 
current,  perhaps  in  different  temperature  ranges. 

2.  Conduction  by  electrons  contributed  by  impurities. 

It  would  then  be  coincidental  that  Cs I  from  Harshaw 
and  Isotopes,,  Inc.  give  similar  results. 

a.  Conduction  by  electrons  produced  at  the  electrodes 

by  photo-  or  the rma 1 - emi s~ l on  due  to  the  hot  specimen. 

This  is  not  a  possible  conduction  mechanism  because 
the  electrons  are  not  free  to  move  through  the  crystal; 
it  might  however,  play  a  role  in  the  polarization 
phenomenon . 
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VII  CONCLUSIONS 


The  alkali  halides,  Csl  and  KI ,  experience  rapid  increases  in  elec¬ 
trical  conductivity  when  subjected  to  shock  waves  of  amplitude  greater 
than  165  kbar  and  185  kbar  respectively.  The  corresponding  shock  tem¬ 
peratures  are  estimated  to  be  2470°K  and  2600  k.  The  results  of  trans¬ 
verse  geometry  shots  with  Csl  (and  KI — see  Report  I)  suggest  that  the 
conductivity  is  dependent  on  the  current  direction  relative  to  the  direc¬ 
tion  of  shock  propagation  (see  Section  IV).  The  discrepancy  between  the 
Csl  data  from  Report  I  and  the  present  data  suggest  also  a  dependence  on 
crystal  orientation  and'or  crystal  purity. 

NaCl  shocked  to  comparable  pressures  has  a  much  lower  temperature 
(due  primarily  to  its  large  specific  heat),  and  its  resistivity  is 

104  O-crn  at  225  kbar.  If  the  specimen  is  preheated  to  480°C  and  shocked 
to  244  kbar,  the  resistivity  is  ~  470  fl-cn.  This  is  the  only  shot  for 
which  we  were  able  to  observe  appreciable  conductivity  in  NaCl. 

We  have  not  yet  determined  the  conductivity  mechanism.  The  data  are 
not  consistent  with  either  of  the  two  models  considered — an  intrinsic  semi¬ 
conductor  with  band  gap  narrowed  by  compression  of  the  lattice,  or  a  simple 
ionic  solid  which  is  melting  at  the  high  temperatures  produced  behind  the 
shock  front. 
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VIII  SUGGESTIONS  FOR  FURTHER  WORK 

1.  Extend  the  measurements  to  specimens  compressed  along  crystal¬ 
lographic  directions  other  than  those  presently  employed. 

2.  Make  additional  conducti"  .  _y  measurements  in  a  transverse 
geometry,  i.e.,  in  a  direction  parallel  to  the  shock  front.  Preliminary 
results  indicate  a  directional  dependence  of  the  conductivity. 

3.  Extend  conductivity  measurements  on  material  in  different  P-V-T 
states,  using  specimen  preheating  and  shock  reflection  techniques. 

4.  Study  crystals  in  which  the  character  and  concentration  of  de¬ 
fects  have  been  altered  by  irradiation,  plastic  deformation,  or  doping. 

5.  Extend  work  on  NaCl  to  higher  temperatures  and  compressions. 

6.  Study  the  polarization  signal  to  determine  its  origin  and  to 
reduce  the  uncertainty  in  interpreting  the  conductivity  records. 


j 
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ESTIMATION  OF  SHOCKED  STATES 
IN  PREHEATED  SPECIMENS 
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APPENDIX  A 


ESTIMATION  OF  SHOCKED  STATES 
IN  PREHEATED  SPECIMENS 

Starting  from  the  Hugoniot  with  initial  state  (F0,T0),  we  will  esti- 
mate  the  Hugoniot  for  initial  state  {V0,T0),  where  the  primed  initial  state 
is  obtained  from  the  unprimed  initial  state  by  heating  at  P  -  0  (Figure  A-l). 


FIG.  A-l  ILLUSTRATING  CALCULATION  OF  HUGONIOT 

OF  PREHEATED  MATERIAL  (1)  FROM  HUGONIOT 
OF  ROOM  TEMPERATURE  MATERIAL  (2) 


If  /3  =  thermal  expansion  coefficient, 

I 

To 

V'0  =  V0  +  f  VfidT  =  V0[l  +  /3{T'0  -  T0 )]  (A-l) 

T  o 

We  will  assume  that  (BP/BF)^  =  constant.  That  is,  in  the  Mie-Griineisen 
form  of  equation  of  state 

r 

P2  ~  Pl  -  7  {E2  ‘  £1>  ’  <A*2) 
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w«‘  wi  l  I  assume  I fV  -  (I’0/P  )  -  constant.  This  isa  drastic:  assumption,  hut 
scorns  to  be  abon  t  a  s  v  a  1  i  d  a  volume  dependence  Tor!' as  that  determined  by  other 
means . 


At  any  (P,  P) 
internal  energy  E 


point  on  a  liugoniot  with  initial  state  ( /*  0 ,  /  0  )  the 
is  given  in  terms  of  the  initial  internal  energy  E ,  by 


E  -  E, 
0 


\  IP  +  P*)(V 


o 


V) 


(A- 3) 


Letting  PQ  =  0,  we  now  apply  Eqs.  A3  and  A- 2  to  determine  the  energy 
at  point  C  by  two  different  paths: 

1 

Path  ADC:  E  =  E0  +-  Pl{VQ  -  V')  +—  (P'  -  Pj)  .  ( A- 4 ) 

2  ^0 


Path  ABC:  E'  =  E'  t  -  P'(V'  -  V' )  =  En  + 

u  2  o  o 


0  i 

C  dT+-  P'(V'  -  P') 


Using  Eq.  A-l, 


E'  =  Fr,+cn(r;-rft)  +  7  P'1(P0-P')  +  p/ir;  -  r0)]  .  (a-5) 


o  ' ‘ o  'o'  2 


and  (A-5)  and  solving  for  P'  yields 

p,fi  *,<i  -  v/v0)  -  vr.l  -  c/r;  -  T.) 

■  ft - r~r~- -  (A_6) 

[2  %(1  -  ‘"''V  -  Vr.J  +  2  -  V 

and  T'  =  r,  +  (£'  -  *,)/C,  =  T,  +  V0(P’  -  P,)/r0C,  .  (A-l) 

Then,  given  the  material  constants  ,  Cp ,  /3,  and  the  Hugoniot 
through  (PQ  =  0,  PQ),  we  can  calculate  the  pressure  P'  and  temperature 
T‘  along  a  Hugoniot  through  (PQ  =  0,  Pg ) •  The  procedure  is  to  assume 
a  P'  ,  calculate  P'  from  Eq.  ( A- 6 ) ,  and  then  T'  from  Eq.  (A-7). 


Equating  ( A- 4 ) 


P' 


To  determine  the  shocked  states  actually  reached  in  a  given  system,  it 
is  necessary  to  map  the  new  P,  V  Hugoniot  into  the  P-u  plane.  The  particle 
velocity  u'  corresponding  to  ( P ' ,  P' )  is  given  by 


u'  =  [P'(P 0  -  P')]'/j 


( A-  8 ) 
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APPENDIX  B 


PRESSURE  DEPENDENCE  OF  IONIC  CONDUCTIVITY 


At  P  *  0,  conduction  in  the  alkali  halides  is  completely  ionic;  the 
mechanism  for  this  is  the  diffusion  of  vacancies.17  The  conductivity 
depends  on  the  density  n  of  carriers  and  on  the  mobility  y.  (the  drift 
velocity  per  unit  electric  field).  We  will  identify  n  with  the  concen¬ 
tration  of  vacancies  of  one  sign,  because  generally  the  mobility  of  one 
ion  is  significantly  greater  than  that  of  the  other. 

For  a  pure  crystal,  there  is  an  equilibrium  value  of  n  (nonzero  for 
T  >  0)  determined  by  the  free  energy  of  formation  A Gj  of  a  vacancy  pair: 


“AG  . / 2k  T 
n  oc  e  ’ 


(&- 1 ) 


The  factor  of  two  arises  because  charge  neutralization  requires  that 
vacancies  be  formed  in  pairs,  and  the  total  number  of  ways  to  choose  a 
pair  of  vacancies  is  the  square  of  the  number  of  ways  to  produce  a  single 
vacancy.  Even  for  the  highest  purity  crystals  available,  n  is  governed  by 
Eq.  B-l  only  for  temperatures  above  perhaps  300°C,  called  the  intrinsic 
region.  Below  this  temperature  is  the  extrinsic  region  where  n  is  deter¬ 
mined  by  certain  impurities.  The  most  important  of  these  impurities  are 
those  with  a  valence  of  +2  (e.g.,  CaCl 2  in  a  NaCl  lattice)  because,  for 
each  divalent  ion  which  takes  a  regular  lattice  position,  there  must  be 
a  positive  ion  vacancy  to  maintain  charge  neutrality. 

The  mobility  is  determined  by  the  jump  frequency  I-1,  which  depends 
on  the  difference  in  free  energies  A G_  of  the  activated  state  (ion  midway 
between  two  lattice  sites)  and  the  ground  state  (ion  at  a  normal  lattice 
site).  The  relation  is  of  the  form 


1  “AC  /k  T 

ll  cc  ~  az  ~  e 

T  T 


(B-2) 
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riu>  exact  form  of  llie  pre -exponent  i  a  1  factor  depends  on  tin-  model  used 
in  the  calculation.  The  conductivity  is  then 


-AC/* T 


(B-3) 


where  A G  =  1/2  A +  A G*.  For  comparison  with  periment,  the  temperature 

dependence  of  the  coefficient  is  often  neglected  and  the  expression 


cr  =  o  e 
00 


-AC/*  T 


(B-3a ) 


is  written. 


The  pressure  dependence  of  cr  can  be  expressed22  by  an  activation 
volume  AF  defined  by  AF  =  (BA G/^P)^.  Differentiating  Eq.  B-3a: 

A*  Jk  I.  cr_\  ai; 


/3  In  cA  /  d  In  CTo>\  AF 

\  V  i  =  l  ^  )T  "kT 


(B-4 ) 


2  o  •  3 

Now  cr^  ccNalv,  where  N  -  number  ions  of  given  sign/cm  cc  F  ,  a  =  lattice 

constant  cc  F1^,  and  v  =  vibration  frequency.  Let  J  oc  8,  the  Debye  tem¬ 
perature,  so  that  v  cc  F  ?  where  y  is  the  Griineisen  parameter  defined  by 
y  -  -( d  In  1 9)/(d  In  F).  Therefore 


cr  ccZVaZ*,  ^  y- 3+2/J-y 


(B-5 ) 


-(’•SC-iH  •(»!)'■. 


(B-6) 


where  B?  is  the  isothermal  bulk  modulus. 

The  ratio  of  the  magnitudes  of  the  two  terms  on  the  right  side  of 
Eq.  B-4  is  then 


(y  1 

B^V 


(B-  7 ) 
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At  P  =  0  and  room  temperature,  y  -  2,  By  -  7  x  10* 'dyne/cm2 , 

AF  =  (AV)fl  =  volume/ion  =  1  *  l(T23cm3.  Then  at  P  =  0,  7  =  300°K, 
r  -  r 0  =  0.03. 

from  the  work  of  Keyes23  we  deduce  that  r  should  never  get  much 
larger.  Keyes  concludes  that  (approximately) 


Bt  Ay  =  sAG 


(B-  8 ) 


Since  s  =  constant  and  y  changes  little,  r  oc  7/AG.  In  integrated  form, 
Eq.  B- 8  becomes 


AG  =  (AG) 


a' 


(B-9) 


which  is  strictly  true  only  at  constant  7.  However,  it  is  consistent  with 
Eq.  B-8  to  assume  that  AG  is  a  function  only  of  volume.  Thus 

r  oc  TV*  (B- 10) 


T  cannot  increase  by  more  than  a  factor  of  ~3  at  P  =  0  without  melting 
occurring.  Increasing  P  enables  T  to  be  increased  but  V  decreases  and, 
since  s  is  in  the  range  of  2  to  4,  r  should  increase  little,  if  any. 

This  argument  is  presented  to  justify  dropping  the  first  term  on  the 
right  side  . of  Eq.  S-4.  Doing  so  and  integrating  gives 


cr  =  cxTe 


•  P&v/kT  [-(AG)0-/>at'] /*7 

=  (7  £ 

co  w 


(B- 11) 


—  _  1  ,p 

where  AE  =  —  J  A VdP  and  subscript  0  indicates  P  =  0. 


Equation  B-ll  provides  an  atomistic  approach  to  the  consideration  of 
the  pressure  dependence  of  & ,  because  the  activation  volume  can  be  inter¬ 
preted  in  terms  of  atomic  volumes  with  fair  success. 22 


Our  problem  is  to  make  a  reasonable  estimate  of  an  upper  limit  for 
ionic  conductivity  at  pressures  and  temperatures  far  exceeding  the  range 
of  available  measurements.  Since  cr  increases  with  T,  we  will  confine  our 
analysis  to  the  region  of  the  melting  curve. 
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Our  first  approach  was  to  use  Eq.  B- 11,  making  the  arbitrary  assump¬ 
tion  that 

V 

AV  =  —  (AV)n  (U-12) 

v  0 

Equation  B-ll  can  be  expressed  in  terms  of  the  volume  along  the  melting 
curve  by  using  the  Simon  expression24  to  relate  P  and  T, 


P  = 


and  the  Murnaghan  expression25  to  relate  P  and  V 


P 


B 


-  BT/b 


(B- 13 ) 


(B-14) 


Both  Eqs.  B-13  and  B-14  are  here  treated  as  empirical  relations.  Combining 
Eqs.  B-ll  and  B-14  gives 


(AG), 


kT 


•  o 


1  + 


B(  AV00 

c1ml 


(B-15) 


where  =  (B/A)l/a,  C2  =  (A/B)~  1. 


In  searching  for  an  improvement  over  the  assumption  of  Eq.  B-12,  the 
relation  (Eq.  B-9)  of  Keyes  was  discovered.  Combining  Eqs.  B-9,  B-ll, 
B-13,  and  B-14  gives 


(B-16) 


^  a*  . 

Equations  6-15  and  B-16  have  been  evaluated  for  sodium.  The  constants 
B  and  6  are  from  fits  to  Bridgman’s  isothermal  data,27  but  the  variation 
from  the  melting  curve  should  not  be  great.  The  constants  used  are  given 
in  Table  B- I  and  the  results  in  Table  B-II.  The  value  of  1.5  for  V 0/T 
corresponds  to  ~100  kbar  in  Na. 
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T.i  l>  1 1-  H- 1 

(  (Asrwrs  rsKI)  IN  KVAIJ'AT'.NU 
K^l  Ml  (AS  H-ir.  AM)  ll-ii,  KOH  SODIUM 


l  a  1,1.  15-  I  I 

CAIjClilATKI)  HKSIJLTS  KOH  SODIIAI 


QUANTITY 

VALUK 

HKKEHKNCE 

A 

11.4  kbar 

28* 

a 

3.56 

28* 

B 

28  kbar 

27b 

b 

3 

27 

T 

mo 

(aK)0 

371°K 

12.2  ctnVmole 

26 

(aG)0 

7800  cal/mole 

26 

s 

2.1 

23 

*  These  constants  were  derived  by  Gilvarry  for 
1 0  kb  data,  but  reproduce  fairly  well  the 
recent  data  of  Newton  et  al.  to  SO  Itbar, 

6  Chosen  to  give  Aaier,  Inst,  Phya,  Hdbk,  value 
for  Bj  (-  bB ). 


vQ/v 

[<AC>0/*  Tm(}] 

Eq.  B-15 

HD 

1 

1 

1 

l.i 

17 

1.4 

1.2 

160 

2.9 

1.5 

30,000 
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The  result  is  that  the  Keyes  expression,  Eq.  B-9,  with  the  experi¬ 
mentally  determined  exponent,  indicates  appreciably  less  variation  of 
conductivity  along  the  melting  curve  than  does  the  arbitrary  assumption 
of  Eq.  B- 12. 

The  possibility  was  explored  that  a  law  of  corresponding  states 
might  apply  to  ionic  conduction.  The  difficulty  here  is  that,  although 
activation  volumes  have  been  determined  for  a  number  of  salids,  only  in 
one  or  two  cases  has  the  pressure  been  raised  high  enough  to  measure  any 
volume  dependence  of  AK,  26,29  Available  data  do  suggest  the  following:2^21 


A H 


T 


m 


constant  for  a  given  material 


(B-17) 


Writing  A G  =  A H  ~  T&S  in  Eq.  B-3a  and  with  T  -  l'm- 

„  ,  ,  -A H/k T 

a  =  e  -  .  (B-18) 

On  the  basis  of  Eq.  B-17,  all  the  volume  dependence  is  in  the  entropy  term, 
which  is  small  6  and  expected  to  vary  slowly. 

A  different  approach  is  to  look  at  the  pressure  dependence  of  dif¬ 
fusion  or  ionic  conductivity  in  the  molten  state  .  Data  have  been  found 
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for  only  two  materials,  Hg  30  and  Ga  .31  In  both  cases,  the  activation  volumes 
for  diffusion  are  of  the  order  of  5  percent  of  atomic  volumes,  much  smaller 
than  expected  for  solid  material.  That  is,  the  pressure  effect  was  found 
to  be  very  small  (ar  is  also  the  case  for  the  temperature  effect)  relative 
to  that  in  solids.  Now  the  conductivities  of  several  molten  salts  have 
been  measured  at  zero  pressure  by  Yaffe  and  Van  Artsdalen20  and  found  to 
be  of  the  order  of  l(Q-cm)-1.  The  temperature  effect  was  very  small — the 
activation  energies  (which  varied  somewhat  with  T)  were  in  the  range  of 
0.1  to  0.2  ev  for  NaCl,  Csl,  and  KI,  one  order  of  magnitude  lower  than  is 
observed  in  the  solid  state.  Intuitively,  we  might  then  expect  that  for 
low  enough  pressures,  the  effect  of  pressure  on  the  conductivity  of  molten 
salts  would  be  small  also.  This  is  suggested  by  the  observations  on  Hg 
and  Ga  and  by  Eq.  B-8  if  it  is  applicable  to  liquids.  This  latter  question 
has  been  considered  by  Keyes24 and,  from  an  examination  of  viscosity  data 
for  isopropyl  alcohol  and  a  silicone  trimer,  he  concludes  that  Eq.  B-8  is 
satisfied  only  below  5  to  10  kbar.  At  higher  pressure  AV  increased  rather 
than  decreased.  This  may  merely  indicate  that  activated-state  theory  is 
not  applicable  to  viscosity  in  liquids. 

In  any  event,  we  expect  the  resistivity  of  a  molten  salt  to  increase 
with  increasing  pressure,  although  the  effect  may  be  small.  Hence,  at 
high  pressure  in  the  vicinity  of  the  melting  curve,  under  equilibrium 
conditions,  we  would  expect  the  alkali  halides  to  have  resistivities 
>  1  Q-cm  even  if  molten,  and  still  greater  if  solid. 
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APPENDIX  C 


APPROXIMATE  CORRECTION  FOR  NONUNIFORM  SHOCK 


The  explosive  systems  used  in  the  present  experiments,  being  finite 
in  size,  do  not  produce  uniform  shock  waves;  in  other  words,  the  pressure 
drops  off  behind  the  shock  front  rather  than  remaining  constant.  The  rate 
of  decay  increases  with  increasing  peak  pressure. 

It  is  only  with  the  recent  development  in  these  laboratories  of  the 
Manganin-wire  shock  pressure  gauge  that  a  measurement  of  pressure  profile 
has  been  possible  Peak  pressures,  as  determined  from  free-surface  velocity 
and/or  impedance  match  measurements,  were  reported  in  Report  I;  in  the 
present  work,  an  attempt  has  been  made  to  correct  approximately  for  the 
pressure  release  behind  the  shock  front.  The  procedure,  chosen  for  sim¬ 
plicity  consistent  with  the  limited  available  data  on  pressure  profiles  and 
the  uncertainty  in  calculating  off-Hugoniot  points,  is  described  below. 

For  each  explosive-driver  plate  assembly,  a  typical  value  of  dP/dt 
(behind  the  front)  was  determined  from  gauge  measurements  made  essentially 
at  the  surface  of  the  driver  (the  gauge  wire  was  actually  displaced  1  mm 
from  the  surface).  Using  this  slope  and  neglecting  any  attenuation  of  the 
peak  pressure  P^,  an  average  specimen  pressure  P'  was  deduced,  corresponding 
to  the  time  of  measurement  of  resistance  from  the  oscilloscope  records.  In 
order  to  estimate  the  volume  and  temperature  corresponding  to  P' ,  the  isen- 
trope  along  which  the  expansion  occurs  was  approximated  by  the  Hugoniot  (see 
Figure  C- 1 ) .  The  volume  corresponding  to  P'  is  then  and  the  corresponding 
temperature  T 2  is  given  by 


*  „  -a(V-V  ) 

r2  =  Tj  e  21 


(C-l) 


where  a  =  y/Y,  y  being  Griineisen’s  ratio,  and  it  is  assumed  that  "a”  is 
constant  ( =  70/FQ ) .  * 


Equation  C-l  expreaaea  the  variation  of  taaiperature  with  voluaie  along  any  iaentrope;  in  the  preaent  ct.ee 
the  iaentrope  through  haa  been  approxiaated  by  the  Hugoniot.  The  leaperature  7^.  correaponding  to 
isentropic  expanaion  froa  Vj  to  Kj*  will  alwaya  exceed  the  teaperature  7*.  on  the  Hugoniot  at  voluae  V 
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•A-  4100- 19 

FIG.  C-i  ESTIMATION  OF  POINT  ON  RELEASE  ISENTROPE 

Now,  as  indicated  in  Figure  C- 1 ,  the  “true”  final  state  is  P'  ,V'  ,T'  , 
and  in  the  following  discussion  we  estimate  this  state  for  comparison  with 
P2  =  P‘  ,V2,T 2  .  The  worst  case  has  been  chosen  —  the  highest  pressure  state 
which  was  produced  by  9404  high  explosive. 

The  problem  is,  given  a  Hugoniot  curve  along  which  P,  V,  and  T  are 
known,  to  determine  V'  and  T'  corresponding  to  a  given  P'  ,  the  primed  point 
lying  on  the  isentrope  through  the  point  P1 , Kj , Tj  on  the  Hugoniot  (see 
Figure  C-l).  Two  equations  in  the  unknowns  V'  and  T'  are 

,  -0(  v‘ -V,  ) 

7  =  i’j e  1  (C-  2) 

P'  -  PZ(V‘)  =  (Bp/Bj) v [T'  ~  T3(V')]  (C- 3) 

We  will  assume 

(BP/BP)  =  aCv  =  constant 


i 


Eliminating  T‘  from  Eqs.  C-2  and  C-3  yields 

P'  ~  De~aV‘  =  P3(V‘)  ~  AT3(V')  ( C- 4 ) 

where 

A  =  aC  and  B  s  AT,e  1 

V  1 

Equation  C-4  was  solved  for  V  by  plotting  the  left  and  right  sides 
of  the  equation  against  V'  as  a  variable  and  finding  the  intersection  of 
the  resulting  curves.  T'  was  then  determined  from  Eq.  C-2.  A  comparison 
of  values  for  the  worst  case  ( i.e .  ,  where  the  corrections  are  largest)  is 
as  follows: 


Peak  values:  P , 

=  340  kbar  , 

*1 

-  0.1290  cm3/g  , 

Tr 

=  5780°K 

“Corrected”  values:  P‘ 

=  298  kbar  , 

=  0. 1325  cm3/g  , 

=  5600 °K 

(T2 

=  4630 °K) 

“True”  values:  P‘ 

=  298  kbar  , 

V 

=  0.1346  cm3/g  , 

r 

=  5500 °K 

Thus  the  volume  correction  is  about  38  percent  too  small  and  the 
temperature  correction  about  36  percent  too  small.  That  is,  the  corrected 
values  of  volume  and  temperature  are  ~1 . 6  percent  too  low  and  1.1  percent 
too  high,  respectively.  The  latter  discrepancy  is  insignificant  with  in¬ 
spect  to  the  absolute  uncertainty  in  calculated  temperatures. 
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APPENDIX  D 


DERIVATION  OF  EQUATION  (9)  FOR  RESISTIVITY 
OF  A  POUR  SEMICONDUCTOR 


For  an  intrinsic  semiconductor  with  equal  concentrations  n  and  equal 
mobilities  of  electrons  and  holes,  the  conductivity  cr  is  given  by 

cr  =  2  nefj.  (D-i) 

where  e  is  the  electronic  charge  and  fj.  the  mobility.  Writing  explicitly 
the  expression  for  n  in  Eq.  (D-l)  yields 

cr  =  le(27TmkT/h2)3/2e~Et/2kT^  ( D-  2 ) 

For  polar  semiconductors32  i t  is  found  that  the  mobility  depends  on 
temperature  as 


p  =  p0  r*  (d-3) 

where  n  usually  lies  between  2  and  3;  let  n  =  2.5.  We  will  evaluate 
by  fitting  Eq.  (D-3)  to  a  typical  value  of  mobility,  viz.,  100  cm2/v-sec 
at  T  =  500°K.32  Hence  pQ  =  5.6  x  108  cm2deg2  •  5 /v- sec .  Then 

cr  =  (B/T)e~E*/2kT  (D-4) 

where 

b  =  4efi0(2irmk/h2)3''2  =  0.85  x  106  “K/D-cm 

Or,  since 

P  =  1  /'cr  ,  P  =  (1.2  x  lO'6)TeEe/2iT  (D-5) 
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APPENDIX  E 


DISCUSSION  OF  PHASE  CHANGES  IN  SODIUM  CHLORIDE 


Evdokimova  and  Vereshchagin32  have  recently  reported  a  crystal  struc¬ 
ture  change  in  NaCl  at  17  kbars  using  high-pressure  X-ray  methods.  They 
report  the  high  pressure  phase  has  a  CsCl  structure  with  a  volume  change 
of  14.9  percent  and  that  no  more  than  5  percent  of  the  substance  transforms 
during  10  to  20  hours  at  '^19  kbar.  They  also  state  that  the  high-pressure 
phase  is  usually  preserved  after  the  release  of  pressure.  These  facts  in¬ 
dicate  an  extremely  sluggish  transition. 

Pistorius34  states  that  Piermarini  and  Weir  were  unsuccessful  in 
corroborating  the  transition,  possibly  due  to  its  sluggishness.  Bassett 
and  Takahashi 35  also  failed  to  corroborate  the  transition.  However, 
Pistorius34 has  succeeded  in  verifying  this  transition  and  in  addition  has 
measured  its  temperature  dependence  to  220°C.  34  He  indicates  that  the 
transition  is  sharper  at  the  higher  temperatures. 

Recently  Larson  at  the  Lawrence  Radiation  Laboratory  at  Livermore36 
has  observed  what  he  believes  to  be  a  multiwave  shock  front  in  single 
crystal  NaCl  using  a  quartz  transducer.  He  is  tentatively  suggesti ng  tha t 
this  multiwave  structure  is  due  to  the  transition  observed  by  Pistorius. 

His  preliminary  data  place  the  transition  at  29  kbars. 

Alder13  has  published  data  indicating  a  break  in  the  Hugoniot  which  he 
suggested  was  a  transformation  from  the  NaCl  to  the  CsCl  type  structure. 

The  pressure  of  the  transition  was  originally  given  as  approximately 
300  kbar.  However,  a  recent  revision  of  the  data36  has  placed  this  trans¬ 
ition  at  about  220  kbar.  The  exact  pressure  seems  to  differ  by  10  or 
20  kbar  depending  on  the  crystal  orientation.36 

This  transition  nas  been  confirmed  qualitatively  for  the  [  1 00 ]  direc¬ 
tion  by  Hauver  and  Melani37who  are  presently  studying  it  as  a  function  of 
crystal  orientation. 
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If  the  preliminary  observation  by  Larson  is  substant  i  ute<l ,  then  two 
transitions  occur  in  NaCl  under  shock.  Larson  suggests  that  the  transition 
at  about  240  kbar  may  be  melting.36 

This  melting  hypothesis  is  not  consistent  with  extrapolation  of  the 
melting  curve  shown  in  Figure  21,  which  indicates  a  melting  temperature  at 
250  kbar  of  about  3000°K;  the  calculated  shock  temperature  at  250  kbar  is 
only  about  1200°K.  However,  the  validity  of  such  a  long  extrapolation  of 
the  static  melting  data  is  highly  questionable;  indeed,  if  the  low-pressure 
transition  is  occurring,  the  extrapolation  is  for  the  wrong  phase. 
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